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CHAPTER I 


INTRODUCTION 

This study analyzes many of the standard stability index 
concepts used today in the prediction of convective weather, 
including severe thunderstorms and tornadoes. These indices are 
computed for the National Aeronautics and Space Administration 
(NASA) Atmospheric Variability Experiment (AVE) mean atmospheric 
soundings (vertical profiles) taken during the AVE-IV project in 
1975 [1]J The profiles have been categorized to correspond to the 
differing weather conditions that occurred by using the manually 
digitized radar data (MDR) taken during this time period as being 
representative of the various weather categories that existed during 
AVE-IV. The AVE procedures are discussed in more detail in 
Chapter III. It is the purpose of this study to present and analyze 
the vertical weather profiles taken during AVE-IV, in terms of 
averaged (mean) weather profiles. 

This investigation also compares several stability indices 
calculated from the AVE-IV mean profiles. The results and per- 
formance of the indices are discussed. 

Also, it is the intent of this study to determine if averaged 
weather profiles taken three hours prior to severe weather occurrence 


Numbers in brackets refer to similarly numbered references 
in the Bibliography. 



have forecast capability when based only on the AVE-IV three-hour 
lag data soundings. 

Finally, the results obtained from this stability index/ 
AVE-IV study are applied to the independent AVE-SESAME-I [2] data 
set to see how well the conclusions based on the AVE-IV analysis 
apply to a set of independent, individual sounding profiles that 
deal with a similar weather system. 


2 



CHAPTER II 


DEFINITIONS 
A. Introduction 

This chapter presents a discussion of the Skew-T, Log-P 
thermodynamic diagram which was used extensively throughout the 
study and in the presentation of results. Secondly, a brief dis- 
cussion defines the stability index. This is concluded by a section 
devoted to describing the different definitions of the stability of 
the atmosphere. 


B. Skew-T, Log-P Diagram 

Throughout this investigation, stability indices, atmospheric 
processes, and atmospheric analyses will involve the use of a suit- 
able meteorological thermodynamic diagram. This will better 
illustrate and describe atmospheric processes, as well as allow 
graphical computations. This section presents a brief but 
instructive discussion of the thermodynamic diagram used in this 
study. 

The Skew-T, Log-P diagram is employed in this study because 
it is most widely used in the United States. In particular, the 
United States Air Force (USAF) Air Weather Service (AWS) uses this 
diagram exclusively. Its popularity is due to the convenience and 
ease of use for most atmospheric computations. The diagram will 
hereafter be referred to as the "Skew-T." The Skew-T contains the 
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same meteorological parametric lines as other thermodynamic diagrams, 
but their arrangement is different. Further discussion of the 
advantages obtained in using the Skew-T diagram is given in 
Reference [3]. 

The name of the diagram indicates how the pressure and 
temperature lines are presented. Constant pressure lines (isobars), 
expressed in millibars (mb), are plotted horizontally on a loga- 
rithmic scale. Also, constant temperature lines (isotherms), in °C, 
are plotted sloping from the lower-left to the upper-right (skewed) 
on an approximate 45 deg angle with respect to the horizontal 
pressure lines. Figure 1 illustrates the isobaric and isothermal 
lines on a Skew-T diagram. 

Also plotted on the Skew-T are slightly curved, dry 
adiabatic lines (°C). They slope from lower-right to upper-left and 
are indicated by two small and one large alternating dashed lines. 
These lines indicate the rate of temperature change encountered when 
a parcel of unsaturated air rises or descends adiabatical ly (without 
gain or loss of heat). 

Saturation adiabats (or moist adiabats, °C) are the large- 
dashed, slightly curved lines also sloping from lower-right to 
upper-left. They begin more vertically at lower levels on the chart 
and become more parallel with the dry adiabatic lines at higher 
levels. Moist adiabatic lines represent the path along which a 
saturated air parcel rises. Figure 1 shows the placement of the dry 
and saturated adiabatic lines on the Skew-T diagram. 


4 










Finally, saturation mixing-ratio (w $ ) lines are shown as 
dashed, slightly curved lines extending from lower-left to upper- 
right in Figure 1. The mixing ratio of an air sample is a function 
of temperature and pressure. It is defined as the ratio of the mass 
of water vapor to the mass of dry air containing the vapor (gm/kg). 
At a given isobaric level, the intersection of the temperature line 
with the w $ line gives the saturation mixing ratio value of the air 
at that temperature and pressure. The dew-point intersection with 
the w $ line gives the actual mixing ratio value (w) of the air. 

To illustrate the use of the Skew-T, an example sounding 
profile of temperature (T) and dew-point temperature (T^) from the 
1,000-mb level upward is plotted on the Figure 2 diagram. Dry 
adiabatic lifting of a surface air parcel is assumed to take place 
in this example. Beginning at the intersection of the T^ and 
1,000-mb pressure line and following the w $ line upward to where it 
intersects the path of the dry adiabat extending upward from the 
surface value of T, introduces an intersection point on the Skew-T, 
called the lifting condensation level (LCL). At this point, 
saturation conditions exist. Traversing vertically from the LCL 
along the saturation adiabat until it intersects the environmental 
sounding of T defines the level of free convection (LFC) location. 
Above the LFC the parcel of air becomes warmer (less dense) than the 
environmental air around it during this period of travel. Above 
this level the parcel will continue to rise at the moist adiabatic 
rate until it becomes cooler than the environment. This, then, 
defines the equilibrium level (EL). As can be seen in Figure 2, 
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Figure 2. 


Skew-T, Log-P diagram, with lifted parcel example [3] 
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regions of negative kinetic energy which work against the vertical 
motion of the cooler air parcel must be overcome. Likewise, positive 
energy areas enhance the parcel's vertical motion. 

A second example is illustrated in Figure 3. In this example 
it is assumed that a parcel of surface air has undergone thermal 
convection produced from solar-ground heating. The parcel rises dry 
adiabatical ly until reaching its convective condensation level (CCL) 
where it becomes saturated. The CCL is the height of the cumuliform 
cloud bases observed in the atmosphere. The CCL is obtained by 
proceeding upward from the surface T^ value (1,000-mb level) along 
the w s line until intersection with the environmental temperature 
sounding occurs. The equilibrium level (EL) is defined in the same 
manner as indicated earlier. 


C. Stability Index 

At least three main factors are determined to be necessary 
for the formation of convective weather: Instability of the atmos- 

phere, sufficient moisture, and a triggering mechanism which lifts 
and sets the atmosphere in motion [4]. Scoggins [4] concluded that 
vertical motion is always required for thunderstorm development, 
regardless of the degree of potential instability. The instability 
of the atmosphere over a location can be calculated by the use of 
upper air data and a stability/instability index computation. With 
the advent of the radiosonde and its routine use in obtaining upper 
air data, stability indices have been developed and used by man 
since the mid-1940's [5], 
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Figure 



3. Skew-T, Log-P diagram, with, heating parcel example [3]. 
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Temperature, pressure, moisture, and winds can be measured 
throughout the upper atmosphere. These data, together with the 
large data processing ability of modern electronic computers, allow 
the researcher to use the data in testing and determining which 
atmospheric parameters vary, and how much, when convective weather 
occurs. This type of parametric study would likely evolve into the 
establishment of a stability index. Generally, stability indices 
take the form of a difference between parameters, such as 
temperature (T), dew-point temperature (T^), potential temper- 
ature (0), mixing ratio (w), pressure (P), altitude (Z), etc., 
measured at two heights or pressure surfaces. The common, available 
pressure levels generally used in index computations are the 1,000-, 
850-, 700- , and 500-mb levels. 

Stability indices act only as an aid in the forecasting of 
convective weather, by alerting the forecaster to areas of the map 
or soundings which should be examined more closely by other methods. 

D. Stability 

Atmospheric instability is usually defined in terms of con- 
ditional instability, latent instability, and potential instability. 
The definitions are not inclusive, however. 

Conditional instability is defined [6] as, "the state of a 
column of air in the atmosphere when its lapse rate of temperature 
is less than the dry-adiabatic lapse rate, but greater than the 
saturation-adiabatic lapse rate. With reference to the vertical 
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displacement of an air parcel, the air will be unstable if saturated, 
and stable if unsaturated." This is illustrated in Figure 4. 

To explain Figure 4, assume parcel theory [3]. When the 
environmental temperature lapse rate (r g ) lies to the left (PQ ) of 
the dry adiabat (T^) through point P, the atmosphere is said to have 
absolute instability within the vertical region between PQ. If an 
air parcel, originally unsaturated, ascends upward along the dry 
adiabat, it will be warmer (at Q 1 ) than the surrounding environment 
(at Q); thereby, the parcel will tend to continue to rise. 

The reverse situation, indicating absolute stability with 
respect to saturation, is true if the environmental lapse rate is 
located to the right (PR) of the saturated adiabat (r s ) PR'. The 
parcel temperature (at R' ) would then be colder than the environment 
temperature (at R), allowing the parcel to sink and be stable. The 
region between the dry and saturated adiabats indicates the region 
of conditional instability. This means that the parcel is stable if 
not saturated, or unstable if saturated. 

The parcel method mentioned here involves simply an 
unsaturated parcel of air which must be forced to ascend vertically 
along a dry adiabat (r^) until saturated at the lifting condensation 
level. It is then forced to ascend along a saturated (or moist) 
adiabat (r s ) from this point upward through the level of free con- 
vection, and thereafter is accelerated along T g by a positive 
buoyancy and need not be forced. Figure 5 shows this process, with 
the arrows indicating the parcel's path. 
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T— ► 

Figure 5. Latent instability example. 
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Latent instability is defined [6] as, "the state of that 
portion of a conditionally unstable air column lying above the level 
of free convection." The negative region (on Figure 5) shown below 
the LFC is the area in which the environment is warmer than the 
parcel. Therefore, if the parcel is initially given an impulse with 
sufficient kinetic energy to carry it through the negative region, 
then above the LFC lies the positive region which signifies the 
latent instability needed to accelerate the parcel, since the parcel 
will now be warmer than the environment. 

Potential instability (or convective instability) is the 

last atmospheric instability category to be considered here. It is 

defined [6] as, "the state of an unsaturated layer or column of air 

in the atmosphere whose wet-bulb potential temperature (9^), or 

equivalent potential temperature (Q^), decreases with altitude. If 

such a column is lifted bodily until completely saturated, it will 

become unstable" (i.e., r , >r ). In this definition one is 

layer s 

considering the stability of a whole layer of air (not a small 
parcel) which is lifted entirely by either frontal activity or flow 
over a mountain. As shown in Figure 6, the bottom of this layer 
(AB) may saturate, via dry/moist adiabatic processes (at A), before 
the top of the layer does (at B). This results in the layer lapse 
rate (between AB) becoming, in time, an unstable layer lapse rate 
(between A 1 B ’ ) . Potential instability (or stability) is strictly a 
"1 ifted-layer"-type of approach to stability. 
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Figure 6 


T — ► 

Layer potential instability. 
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CHAPTER III 


ATMOSPHERIC VARIABILITY EXPERIMENT IV 
DATA DESCRIPTION 

A. Introduction 

Presented in this chapter is the description of the Atmos- 
pheric Variability Experiment IV (AVE-IV). This includes background 
information for the experiment, the synoptic situation present, the 
data soundings obtained, the corresponding available radar data, and 
the data reduction technique used. Finally, averaged AVE-IV profiles 
pertaining to pre-storm and storm environments are presented for 
different severities of radar-measured weather conditions. 

B. AVE-IV Data Background 

The NASA AVE-IV project [1] took place between 0000 GMT, 

April 24 and 1200 GMT, April 25, 1975. Forty-two AVE network 
rawinsonde stations participated in this 1.5-day mesoscale experi- 
ment in which atmospheric soundings, from the surface to 25 mb, were 
taken at each site every three hours (with some exceptions). 

Releases were taken nine times at most sites: April 24 at 0000, 

0600, 1200, 1500, 1800, and 2100 GMT, and on April 25 at 0000, 0600, 
and 1200 GMT. Figure 7 shows a map of rawinsonde stations, east of 
the Rocky Mountains, that participated in the AVE-IV experiment. 

Table 1 lists each station. Because of the small temporal and 
spatial resolution of these sounding data, it is believed that 
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Figure 7. Location of rawinsonde stations for AVE-IV [11] 
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Table 1. Rawinsonde Stations Participating in AVE-IV Experiment 


Station 

i Number 

Location 

208 

(CHS) 

Charleston, South Carolina 

211 

(TPA) 

Tampa, Florida 

213 

(AYS) 

Waycross, Georgia 

o 

CM 

C\J 

(VPS) 

Apalachicola, Florida 

226 

(CEN) 

Centerville, Alabama 

232 

(BVE) 

Boothville, Louisiana 

235 

(JAN) 

Jackson, Mississippi 

240 

( LCH ) 

Lake Charles, Louisiana 

248 

(SHV) 

Shreveport, Louisiana 

255 

(VCT) 

Victoria, Texas 

260 

(SEP) 

Stephenvi 1 le, Texas 

261 

(DRT) 

Del Rio, Texas 

265 

(MAP) 

Midland, Texas 

304 

(HAT) 

Hatteras, North Carolina 

311 

(AHN) 

Athens, Georgia 

317 

(GSO) 

Greensboro, North Carolina 

327 

(BNA) 

Nashville, Tennessee 

340 

(LIT) 

Little Rock, Arkansas 

349 

(UMN) 

Monette, Missouri 

363 

(AMA) 

Amarillo, Texas 

402 

(WAL) 

Wallops Island, Virginia 

405 

(IAD) 

Sterling, Virginia (Dulles Airport) 

425 

( HTS ) 

Huntington, West Virginia 

429 

(DAY) 

Dayton, Ohio 

433 

(SL0) 

Salem, Illinois 

451 

(DDC) 

Dodge City, Kansas 
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Table 1. (continued) 


Station Number Location 


456 

(TOP) 

Topeka, Kansas 

486 

(JFK) 

Fort Totten, New York (Kennedy Airport) 

518 

(ALB) 

Albany, New York 

520 

(PIT) 

Pittsburg, Pennsylvania 

528 

(BUF) 

Buffalo, New York 

532 

(PIA) 

Peoria, Illinois 

553 

(OMA) 

Omaha, Nebraska 

562 

(LBF) 

North Platte, Nebraska 

606 

(PWM) 

Portland, Maine 

637 

(FNT) 

Flint, Michigan 

645 

(GRB) 

Green Bay, Wisconsin 

654 

(HUR) 

Huron, South Dakota 

655 

(STC) 

St. Cloud, Minnesota 

662 

(RAP) 

Rapid City, South Dakota 

11001 

(MFS) 

Marshall Space Flight Center, Alabama 

22002 

(FSI) 

Fort Sill, Oklahoma 
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smaller meteorological scale (mesoscale) can be studied in terms of 
the variability of atmospheric parameters, than have been studied in 
the past, in application to stability analyses. Normally, across 
the United States, rawinsonde releases take place with a 12-hour 
separation and over a significantly wider spatial network of 
stations. The data reduction and processing procedures, together 
with further project information and the data itself (with 25-mb 
spacing), are presented in Reference [1], 

C. Synoptic Situation 

The surface synoptic weather map for the beginning of the 
AVE-IV experiment (0000 GMT, April 24, 1975) is presented in Figure 8. 
The general weather situation throughout the AVE-IV experiment con- 
sisted of a cold polar air mass moving slowly across the northern 
United States with warm, moist air from the Gulf of Mexico flowing 
over the southern and eastern states. This movement was due to 
circulation around a high-pressure cell located off the coast of the 
Carolinas. At the start of the experiment, these two differing air 
masses were separated by a pseudo-stationary front extending from a 
low-pressure cell over lower Michigan into a secondary low located 
over Kansas. From there, the front trailed into west Texas, as 
shown in Figure 8. Throughout the AVE-IV period, the primary low 
moved into the Gulf of St. Lawrence, while the secondary low had 
moved into Kentucky by the end of the experiment. 

The upper atmospheric flow pattern remained basically zonal 
throughout the experiment, with the exception of two short wave 
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Figure 8. Surface map for 0000 GMT, April 24, 1975 [9]. 
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passages which moved through the network. This wave activity 
resulted in the formation of two squall lines which produced severe 
we at her . 

The first short wave disturbance was already located in the 
Midwest at the beginning of AVE-IV and produced the squall line, 
from Kansas through Illinois, as shown in Figure 8 (at 0000 GMT, 

April 24, 1975). The squall line then moved easterly, ahead of the 
front, and produced maximum thunderstorm activity between 0300 and 
0600 GMT. All thunderstorm, hail, and tornadic activity produced by 
this system had ended by 0000 GMT on April 25, 1975. 

The second short wave passage produced the squall line situ- 
ated through Oklahoma at 0000 GMT, April 25, 1975, as shown in 
Figure 9. Most of the tornadic and severe weather throughout AVE-IV 
resulted from this second squall line as it moved eastward. This 
second squall line formed initially sometime after 2100 GMT on 
April 24, stretching from Missouri into Texas. Storms and convective 
development continued until 0600 GMT, April 25, 1975, when maximum 
squall development occurred, producing large hail, strong winds, and 
tornadoes. This activity included the Neosho, Missouri, tornado at 
0040 GMT, April 25. The line was moving eastward and was still 
strong by the end of the experiment, although the thunderstorm 
activity had lessened. The final surface weather map of AVE-IV for 
1200 GMT, April 25, 1975, is shown in Figure 10. 

The AVE-IV data collection and analyses have been carried 
out by several investigators. Complete AVE-IV information and 
analyses can be found in References [1] and [7 through 22]. 
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Figure 10. Surface map for 1200 GMT, April 25, 1975 [1]. 
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D. Manually Digitized Radar Data 


In order to correlate the stability analyses with the radar 
measurements of precipitation which developed during AVE-IV, the 
manually digitized radar (MDR) data, from the National Oceanic and 
Atmospheric Administration (NOAA) Techniques Development Laboratory, 
were used. These data had been obtained and correlated previously 
for other AVE-IV investigations before use in the present study. 

The MDR grid network of squares (83 km on a side) is shown in 
Figure 11. Areal coverage and echo intensity of rainfall within 
each square, for every hour, determined the MDR code (from 0 to 9) 
assigned to each square. The code used is described by Foster and 
Reap [23] and is given in Table 2. Radar data for each square were 
then compiled, with the maximum hourly radar intensity value over a 
three-hour period being used. The MDR time period was centered on 
each of the nine AVE-IV rawinsonde observation times in order to 
compare the two sets directly. However, for this study, instead of 
using all nine categories of MDR precipitation codes, only four 
composite categories of MDR precipitation intensity/coverage classi- 
fications were used. These MDR definitions were taken from Reap [24] 
and Wilson [17] and are presented in Table 3. 

E. Rawinsonde Data Grid 

The AVE-IV 25-mb spaced rawinsonde profile data were 
interpolated for each of the nine time periods using an 18x18 grid, 
with 160-km spacing between grid points, as shown in Figure 12. 
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Figure 11. Manually digitized radar (MDR) grid network [8]. 


Table 2. Manually Digitized Radar Data Code [8] 


Intensity, Coverage, 

Maximum Percent of VIP Maximum 
Code Observed Coverage Rainfall Intensity 

No. V I P a Values In Box Rate (in h~^) Category 


0 No Echoes 


1 

1 

Any VIP1 


<0.1 

Weak 

2 

2 

<50% of VIP2 

0.1 

to 0.5 

Moderate 

3 

2 

>50% of VIP2 

0.5 

to 1.0 

Moderate 

4 

3 

<50% of VIP3 

1.0 

to 2.0 

Strong 

5 

3 

>50% of VIP3 

1.0 

to 2.0 

Strong 

6 

4 

<50% of VIP3 
and 4 

1.0 

to 2.0 

Very Strong 

7 

4 

>50% of VIP3 
and 4 

1.0 

to 2.0 

Very Strong 

8 

5 or 6 

<50% of VIP3, 
4, 5, and 6 


>2.0 

Intense or 
Extreme 

9 

5 or 6 

>50% of VIP3, 
4, 5, and 6 


>2.0 

Intense or 
Extreme 


a Video Integrator Processor (intensity of returned 
radar signal, gated). 
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Table 3. MDR Categories Used in the Present AVE-IV Analyses 


Category 

A 

B 

C 


MD R Val ue 
0 

>0 

> 3 

> 7 


Convective Activity 

No precipitation 

All precipitation 

All thunderstorm activity 

All severe thunderstorm activity 


D 


4 4 4 
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Figure 12. The 18x18 grid used for numerical computations [8] 







This produced a workable field of measured data at all grid points. 
According to Barr et al . [25], this grid spacing produces the maximum 
resolution possible given a sample of randomly spaced rawinsonde 
stations. More detail concerning numerical computation criteria in 
using the AVE grid is given in the 1976 report of Wilson and 
Scoggins [26]. 

F. AVE-IV Averaged Profiles 

A number of stability analyses have been carried out for 
storm and severe storm environments where stability indices were 
calculated. Most analyses involved the computation of only one or 
maybe two indices. However, there have been only a few studies in 
which a number of stability indices have been computed with identical 
data and compared. Some of these studies are reported in References 
[26 through 31]. Most of these studies involved comparisons of the 
different station indices computed throughout convective development 
of a moving storm system. Also, four reports on stability analyses 
during the AVE-IV project have been published [7, 9, 11, 18]. 

However, as indicated in the Introduction, a different type 
of atmospheric stability analysis will be examined here; that is, 
one involving arithmetically averaged soundings which relate to 
different AVE-IV weather criteria, ranging from no weather to severe 
weather. It was suggested that if mean atmospheric profiles repre- 
senting a certain convective atmospheric environment were compared 
with soundings representing, say, "a more severe environment," then 
an examination of all parameter profile averages might indicate a 
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structural trend within these profiles that would be related directly 
to the degree of convection just prior to or during occurrence of 
severe weather. A forecast tool might result from these trends if 
examined. This average profile study may uncover something unique 
when applied to a convective situation, not observable from an indi- 
vidual station's vertical sounding. It is suggested that atmospheric 
stability through a stability index procedure is one way to do an 
analysis on averaged profile soundings. This averaged-prof i le method 
is not new. Wilson and Scoggins [20] in 1978 presented a quick-look 
at AVE-IV average sounding analyses involving temperature, dew point, 
and vector wind, along with a few calculated parameters. The present 
study extends the work of Wilson and Scoggins [20] in terms of a 
detailed study of just the thermodynamic stability of the AVE-IV 
atmosphere. 

The AVE-IV profile data, related to a grid, can now be linked 
to the MDR grid data. This computational linking had previously 
been done by objective techniques [21] for use in other AVE-IV 
studies. This resulted in producing the six averaged (mean) vertical 
profiles of temperature, dew-point temperature, mixing ratio, zonal 
(east-west) and meridional (north-south) wind speed, and pressure 
level height for the 17 pressure levels of data from 900 to 100 mb, 
with 50-mb spacing for the nine AVE-IV time periods. The procedure 
to obtain these average profiles versus weather category is described 
by example in the following paragraph. 

As an example, consider the most severe thunderstorm cases 
(MDR>7). The following procedure was used to create average 
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soundings for the six measured parameters described in the preceding 
paragraph. All values of the parameters at grid points within ~80 km 
of a three-hour composite MDR value >7 were averaged for the total 
AVE-IV time period. This procedure was also carried out for the 
other three MDR categories defined in Table 3, page 27. The 
resulting four tables of averaged (mean) profiles are presented in 
Tables 4 through 7. 

G. AVE-IV Average Lag Profiles 

As a final task using the AVE-IV data, average lag profiles 
were computed. Lag here is defined as the time difference between 
the sounding and the occurrence of severe weather, three hours later. 
Examination of lag profiles promises a certain forecast capability 
through determination of the average environment three hours prior 
to severe weather occurrence. 

The three-hour composite MDR data had previously been cate- 
gorized according to the four weather types given in Table 3. 

To create the average lag profile, all soundings three hours^ prior 
to the occurrence of each MDR convective category were extracted 
from the data set for each parameter. These were then averaged to 
obtain the average lag profile for the four MDR cases. The results 
are presented in Tables 8 through 11. 


^If three-hour sounding separation was not available, the 
sounding taken six hours prior was used. 
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Table 4. 


AVE-IV Average Profile Conditions for MDR = 0 (No Precipitation) 


Press. 

mb 

Ht. 

m 

Temp. 

°C 

Dew Pt. 
°C 

U Wind 
m/s 

V Wind 
m/s 

Mix. Ratio 
gm/kg 

Wind Sp. 
m/s 

Wind Dir. 
Deg. 

900 

1,010 

13.8 

8.0 

2.0 

4.3 

8.0 

4.7 

205 

850 

1,490 

12.4 

2.9 

4.0 

4.0 

6.3 

5.7 

225 

800 

2,000 

10.1 

-2.6 

5.7 

3.4 

4.7 

6.6 

239 

750 

2,530 

7.6 

-10.5 

7.4 

2.8 

3.1 

7.9 

249 

700 

3,100 

4.4 

-15.5 

8.9 

2.5 

2.2 

9.2 

254 

650 

3,700 

0.5 

-17.9 

10.4 

2.4 

1.9 

10.7 

257 

600 

4,330 

-4.1 

-22.6 

11.9 

2.1 

1.5 

12.1 

260 

550 

5,010 

-9.0 

-26.8 

13.6 

2.0 

1.2 

13.8 

262 

500 

5,740 

-14.2 

-31.8 

15.4 

2.5 

0.8 

15.6 

261 

450 

6,530 

-19.8 

-39.3 

17.4 

3.0 

0.5 

17.7 

260 

400 

7,400 

-26.4 

-44.4 

19.7 

3.3 

0.3 

20.0 

260 

350 

8,350 

-33.7 

-50.4 

22.8 

3.9 

0.2 

23.1 

260 

300 

9,410 

-42.0 

— 

26.2 

4.9 

— 

26.7 

259 

250 

10,600 

-51.1 

— 

30.1 

5.7 

— 

30.6 

259 

200 

12,000 

-59.5 

— 

32.2 

5.1 

— 

32.6 

261 

150 

13,800 

-59.5 

— 

28.2 

5.0 

— 

28.6 

260 

100 

16,400 

-62.2 

— 

19.1 

4.7 

— 

19.7 

256 


Note: Number of soundings = 1,053. 


Table 5. 


AVE-IV Average Profile Conditions for MDR > 0 (All Precipitation) 


Press. 

mb 

Ht. 

m 

Temp. 

°C 

Dew Pt. 
°C 

U Wind 
m/s 

V Wind 
m/s 

Mix. Ratio 
gm/kg 

Wind Sp. 
m/s 

Wind Dir. 
Deg. 

900 

998 

13.7 

9.8 

5.1 

7.6 

8.8 

9.2 

214 

850 

1,480 

11.3 

6.4 

8.1 

7.2 

7.6 

10.8 

228 

800 

1,980 

8.8 

3.1 

10.2 

6.4 

6.4 

12.0 

238 

750 

2,520 

5.9 

-2.4 

11.6 

6.0 

5.0 

13.1 

243 

700 

3,080 

2.6 

-8.0 

12.8 

6.2 

3.6 

14.2 

244 

650 

3,670 

-1.1 

-11.8 

14.4 

6.6 

2.9 

15.8 

245 

600 

4,310 

-5.3 

-17.0 

16.1 

6.8 

2.2 

17.5 

247 

550 

4,990 

-9.9 

-22.2 

18.0 

6.7 

1.7 

19.2 

250 

500 

5,710 

-14.7 

-28.1 

19.7 

6.9 

1.2 

20.9 

251 

450 

6,500 

-20.2 

-34.9 

20.9 

7.4 

0.7 

22.2 

251 

400 

7,370 

-26.5 

-40.9 

22.6 

8.0 

0.4 

24.0 

251 

350 

8,320 

-33.8 

-47.2 

24.9 

8.5 

0.2 

26.3 

251 

300 

9,380 

-42.1 

— 

27.5 

9.6 

— 

29.1 

251 

250 

10,600 

-51.5 

— 

30.5 

10.5 

— 

32.3 

251 

200 

12,000 

-60.6 

— 

32.5 

9.2 

— 

33.8 

254 

150 

13,800 

-59.5 

— 

28.9 

6.7 

— 

29.7 

257 

100 

16,300 

-60.2 

— 

20.2 

4.5 

— 

20.7 

257 


Note: Number of soundings = 567. 


Table 6. 


AVE-IV Average Profile Conditions for MDR > 3 (All Thunderstorms) 


Press. 

mb 

Ht. 

m 

Temp. 

°C 

Dew Pt. 
°C 

U Wind 
m/s 

V Wind 
m/s 

Mix. Ratio 
gm/kg 

Wind Sp. 
m/s 

Wind Dir. 
Deg. 

900 

987 

15.9 

11.2 

5.4 

6.7 

9.6 

8.6 

219 

850 

1,470 

13.3 

7.7 

8.6 

6.6 

8.2 

10.8 

233 

800 

1,980 

10.5 

4.4 

10.8 

6.0 

7.1 

12.4 

241 

750 

2,520 

7.3 

-0.1 

12.3 

5.7 

5.7 

13.6 

245 

700 

3,080 

3.9 

- 7.5 

14.0 

6.2 

3.8 

15.3 

246 

650 

3,680 

-0.1 

-12.1 

15.8 

7.4 

2.9 

17.5 

245 

600 

4,310 

-4.5 

-17.1 

17.9 

8.1 

2.1 

19.7 

246 

550 

4,990 

-9.4 

-22.5 

19.9 

8.5 

1.6 

21.6 

247 

500 

5,720 

-14.3 

-29.1 

21.4 

8.8 

1.1 

23.1 

248 

450 

6,510 

-19.7 

-35.5 

22.5 

9.3 

0.7 

24.4 

248 

400 

7,380 

-26.1 

-40.6 

24.0 

10.0 

0.4 

26.0 

247 

350 

8,330 

-33.5 

-46.6 

25.8 

10.4 

0.3 

27.8 

248 

300 

9,390 

-41.7 

— 

27.8 

11.4 

— 

30.1 

248 

250 

10,600 

-51.1 

— 

30.0 

12.0 

— 

32.3 

248 

200 

12,000 

-60.5 

— 

32.1 

10.9 

— 

33.9 

251 

150 

13,800 

-60.0 

— 

29.3 

8.6 

— 

30.5 

254 

100 

16,300 

-60.7 

— 

20.7 

6.7 

— 

21.8 

252 


Note: Number of soundings = 189. 


Table 7. 


AVE-IV Average Profile Conditions for MDR > 7 (Severe Thunderstorms) 


Press. 

mb 

Ht. 

m 

Temp. 

°C 

Dew Pt. 
°C 

U Wind 
m/s 

V Wind 
m/s 

Mix. Ratio 
gm/kg 

Wind Sp. 
m/s 

Wind Dir. 
Deg. 

900 

978 

18.8 

11.6 

4.3 

6.2 

9.9 

7.6 

215 

850 

1,470 

15.6 

8.2 

7.5 

6.4 

8.5 

9.9 

230 

800 

1,980 

12.5 

5.0 

8.9 

5.6 

7.3 

10.5 

238 

750 

2,520 

8.6 

0.7 

10.2 

4.6 

5.9 

11.2 

246 

700 

3,090 

4.8 

-7.8 

12.1 

5.4 

3.6 

13.3 

246 

650 

3,680 

0.5 

-11.8 

14.4 

7.0 

2.8 

16.0 

244 

600 

4,320 

-4.4 

-17.2 

17.4 

8.0 

2.0 

19.2 

245 

550 

5,000 

-9.4 

-23.5 

19.5 

8.9 

1.4 

21.4 

245 

500 

5,730 

-13.9 

-29.9 

20.5 

9.5 

1.0 

22.6 

245 

450 

6,520 

-19.4 

-35.9 

22.3 

10.0 

0.6 

24.4 

246 

400 

7,390 

-25.9 

-40.4 

24.3 

10.8 

0.4 

26.6 

246 

350 

8,340 

-33.4 

-45.9 

26.7 

10.9 

0.3 

28.8 

248 

300 

9,410 

-41.5 

— 

29.4 

11.8 

— 

31.7 

248 

250 

10,600 

-51.0 

— 

32.8 

12.6 

— 

35.1 

249 

200 

12,000 

-60.2 

— 

35.3 

12.0 

— 

37.3 

251 

150 

13,800 

-60.7 

— 

31.1 

10.1 

— 

32.7 

252 

100 

16,300 

-61.0 

— 

21.8 

7.8 



23.2 

250 


Note: Number of soundings = 66. 


Table 8. 


AVE-IV Average Lag Profile Conditions for MDR = 0 (No Precipitation) 


Press. 

mb 

Ht. 

m 

Temp. 

°C 

Dew Pt. 
°C 

U Wind 
m/s 

V Wind 
m/s 

Mix. Ratio 
gm/kg 

Wind Sp. 
m/s 

Wind Dir. 
Deg. 

900 

1,010 

14.3 

8.2 

2.1 

4.4 

8.1 

4.9 

206 

850 

1,490 

12.5 

3.6 

4.2 

4.1 

6.5 

5.9 

226 

800 

2,000 

10.1 

-1.7 

6.0 

3.5 

5.0 

7.0 

240 

750 

2,530 

7.5 

-9.0 

7.7 

2.9 

3.4 

8.2 

249 

700 

3,100 

4.3 

-14.4 

9.2 

2.7 

2.3 

9.6 

254 

650 

3,700 

0.4 

-16.7 

10.7 

2.7 

2.0 

11.0 

256 

600 

4,330 

-4.1 

-21.3 

12.4 

2.4 

1.6 

12.6 

259 

550 

5,010 

-9.0 

-25.3 

14.1 

2.2 

1.3 

14.3 

261 

500 

5,740 

-14.1 

-30.4 

15.9 

2.6 

0.9 

16.1 

261 

450 

6,540 

-19.7 

-37.7 

17.7 

3.2 

0.5 

18.0 

260 

400 

7,400 

-26.2 

-43.3 

19.9 

3.7 

0.3 

20.2 

259 

350 

8,350 

-33.5 

-49.5 

22.7 

4.5 

0.2 

23.1 

259 

300 

9,410 

-41.9 

— 

26.0 

5.6 

— 

26.6 

258 

250 

10,600 

-51.1 

— 

29.7 

6.6 

— 

30.4 

257 

200 

12,000 

-59.4 

— 

32.0 

6.1 

— 

32.6 

259 

150 

13,800 

-59.3 

— 

28.3 

5.5 

— 

28.8 

259 

100 

16,400 

-62.2 

— 

19.4 

4.7 

— 

20.0 

256 


Note: Number of soundings = 956. 


Table 9. AVE-IV Average Lag Profile 


Press. 

mb 

Ht. 

m 

Temp. 

°C 

Dew Pt. 
°C 

U Wind 
m/s 

900 

999 

13.5 

9.9 

5.3 

850 

1,480 

11.3 

6.1 

8.3 

800 

1,980 

8.8 

2.7 

10.4 

750 

2,520 

6.1 

-3.0 

11.7 

700 

3,080 

2.7 

-8.4 

13.1 

650 

3,670 

-1.0 

-12.5 

14.8 

600 

4,310 

-5.2 

-18.1 

16.4 

550 

4,990 

-10.0 

-22.5 

18.2 

500 

5,710 

-14.8 

-28.2 

20.1 

450 

6,500 

-20.2 

-35.5 

21.4 

400 

7,370 

-26.6 

-41.3 

22.9 

350 

8,320 

-33.9 

-48.0 

25.4 

300 

9,380 

-42.3 

— 

28.0 

250 

10,600 

-51.8 

— 

31.1 

200 

12,000 

-60.6 

— 

32.7 

150 

13,800 

-59.4 

— 

29.4 

100 

16,300 

-60.0 

— 

20.4 


Note: Number of soundings = 484. 


for MDR >0 (All Precipitation) 


V Wind 
m/s 

Mix. Ratio 
gm/kg 

Wind Sp. 
m/s 

Wind 1 
Deg 

8.5 

8.8 

10.0 

212 

7.9 

7.6 

11.5 

226 

6.9 

6.3 

12.5 

236 

6.4 

4.9 

13.3 

241 

6.4 

3.6 

14.6 

244 

6.7 

2.8 

16.3 

246 

6.9 

2.1 

17.8 

247 

6.8 

1.6 

19.4 

250 

7.0 

1.1 

21.3 

251 

7.6 

0.7 

22.7 

250 

7.8 

0.4 

24.2 

251 

8.1 

0.2 

26.7 

252 

9.0 

— 

29.4 

252 

9.6 

— 

32.6 

253 

8.4 

— 

33.8 

256 

6.9 

— 

30.2 

257 

4.7 

— 

20.9 

257 


Table 10. 


AVE-IV Average Lag Profile Conditions for MDR >3 (All Thunderstorms) 


Press. 

mb 

Ht. 

m 

Temp. 

°C 

Dew Pt. 
°C 

U Wind 
m/s 

V Wind 
m/s 

Mix. Ratio 
gm/kg 

Wind Sp. 
m/s 

Wind Dir. 
Deg. 

900 

994 

15.9 

11.5 

5.7 

9.2 

9.7 

10.8 

212 

850 

1,480 

13.3 

8.0 

8.9 

8.5 

8.4 

12.3 

226 

800 

1,990 

10.5 

4.6 

11.3 

7.3 

7.1 

13.5 

237 

750 

2,520 

7.7 

-1.1 

12.9 

6.8 

5.4 

14.6 

242 

700 

3,090 

4.3 

-8.8 

14.5 

7.1 

3.5 

16.1 

244 

650 

3,690 

0.5 

-14.6 

16.6 

7.9 

2.6 

18.4 

245 

600 

4,330 

-4.0 

-19.6 

18.4 

8.6 

1.8 

20.3 

245 

550 

5,010 

-9.0 

-23.9 

20.3 

8.6 

1.4 

22.1 

247 

500 

5,740 

-14.0 

-30.1 

22.4 

8.8 

0.9 

24.1 

249 

450 

6,530 

-19.5 

- 37.6 

23.4 

9.4 

0.5 

25.2 

248 

400 

7,390 

-26.0 

-42.8 

24.5 

9.4 

0.3 

26.2 

249 

350 

8,350 

-33.3 

-49.4 

26.5 

9.3 

0.2 

28.1 

251 

300 

9,410 

-41.8 

— 

28.0 

10.3 

— 

29.8 

250 

250 

10,600 

-51.1 

— 

30.3 

10.4 

— 

32.0 

251 

200 

12,000 

-60.2 

— 

32.4 

9.7 

— 

33.8 

253 

150 

13,800 

-59.9 

— 

30.3 

8.8 

— 

31.6 

254 

100 

16,400 

-60.9 

— 

20.9 

6.2 

— 

21.8 

253 


Number of soundings = 164. 


Table 11. AVE-IV Average Lag Profile Conditions for MDR > 7 (Severe Thunderstorms) 


Press. 

mb 

Ht. 

m 

Temp. 

°C 

Dew Pt. 
°C 

U Wind 
m/s 

V Wind 
m/s 

Mix. Ratio 
gm/kg 

Wind Sp. 
m/s 

Wind Dir. 
Deg. 

900 

992 

18.5 

13.1 

4.6 

8.9 

10.7 

10.0 

207 

850 

1,480 

15.8 

8.9 

7.8 

8.2 

8.9 

11.3 

224 

800 

1,990 

12.7 

4.7 

10.1 

6.3 

7.2 

11.9 

238 

750 

2,530 

9.4 

-0.6 

11.3 

5.4 

5.4 

12.5 

244 

700 

3,100 

5.9 

- 7.7 

12.9 

6.2 

3.5 

14.3 

244 

650 

3,700 

1.8 

-12.6 

15.1 

7.7 

2.5 

17.0 

243 

600 

4,340 

-3.3 

-18.4 

17.5 

9.1 

1.8 

19.7 

243 

550 

5,030 

-8.6 

-22.8 

20.0 

9.3 

1.3 

22.1 

245 

500 

5,760 

-13.6 

-28.0 

23.1 

9.7 

1.0 

25.1 

247 

450 

6,550 

-19.1 

-36.8 

25.1 

10.0 

0.5 

27.0 

248 

400 

7,420 

-25.7 

-41.2 

26.8 

10.0 

0.4 

28.6 

250 

350 

8,370 

-33.1 

-46.1 

29.5 

9.8 

0.2 

31.1 

252 

300 

9,430 

-41.6 

— 

32.4 

10.0 

— 

33.9 

253 

250 

10,600 

-50.7 

— 

35.3 

10.1 

— 

36.7 

254 

200 

12,100 

-59.9 

— 

37.0 

10.9 

— 

38.6 

254 

150 

13,900 

-60.7 

— 

32.6 

10.2 

— 

34.2 

253 

100 

16,400 

-61.4 

— 

22.5 

6.9 

— 

23.5 

253 


Note: Number of soundings = 51. 


CHAPTER IV 


INDICES USED IN STUDY 
A. Introduction 

This chapter presents the criteria used in the selection of 
stability indices that were chosen for analyses in the present study. 
The indices are then presented, with a detailed description given 
for each. 


B. Index Selection Criteria 

The stability indices used in this study were selected to 
utilize the available AVE-IV data described in Chapter III. Indices 
were chosen based on ease of computation. Computations involving 
differences, additions, multiplications, and divisions among the 
available atmospheric parameters at or between vertical pressure 
levels were, in general, selected. 

Mean profile data for AVE-IV does not extend below the 900-mb 
level and therefore, all atmospheric stability indices which use the 
surface or data levels up to 900 mb were eliminated from this study. 
Indices which require complex computation with the available data 
were also eliminated (i.e., indices which require forecasted tempera- 
ture or moisture parameters at the surface or aloft). Finally, 
since the computer was used in computing index values for this study, 
most indices involving a thermodynamic diagram computation were not 
used. 
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C. Indices Chosen 


Fourteen atmospheric stability indices were chosen for 
testing with the AVE-IV mean profile data. They are listed in 
Table 12 and are described in detail in Section D. 

D. Definition of Indices 

This section defines and gives historical information con- 
cerning each stability index used in the study. In order for the 
reader to follow various thermodynamic procedures involved in atmos- 
pheric processes used in the index computation, a simplified Skew-T 
diagram (as described in Chapter II) is given whenever possible to 
help describe and visualize the steps taken during the index compu- 
tation. 

Severe Weather Threat Index 

The Severe Weather Threat (SWEAT) index was developed by the 
United States Air Force Global Weather Central (AFGWC) and presented 
in 1970 [32, 33, 34] for use in forecasting potentially critical 
convective weather (i.e., severe thunderstorms and tornadoes). It 
is a computer-prepared index based on weighted, empirical parameters 
at the 850- and 500-mb levels. The Air Force has revised the SWEAT 
index twice thus far, and all revisions to date will be presented in 
this section. 

The initial SWEAT index (SWEAT -j ) from Reference [32] was 
derived subjectively from a study of 328 severe storm vertical 
soundings and is written as: 


41 



Table 12. Stability Indices 


Index Name 

Symbol 

SWEAT Index 

SWEAT 

Vertical Totals Index 

VTI 

Cross Totals Index 

CTI 

Total Totals Index 

TTI 

Theta E 

0 E 

Showalter Index 

SI 

Rackcliff Index 

RI 

Jefferson Index 

JI 

Modified Jefferson Index 

MJI 

Boyden Index 

BI 

Bradbury Potential Stability Index 

BPI 

K-Index 

KI 

Energy Index 

El 

Modified Martin Index 

MI 
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SWEAT.] 12 T d85Q 


(1) 


+ 20 < TTI - 49 > + 2W 850 + W 500 * 

where, 

T d850 = dew “P oint temperature (°C) at the 850-mb level 
(positive T dg50 values only are used; 

if T d850 < °> then set T d850 = 0) • 

TTI = total-totals Index (°C) 

TTI ' < T + T d>850 ' 2 T 500 
(if TTI <49, set TTI = 49; the second term 

then drops out of Eq. (1)) , 

Wg5Q = 850-mb wind speed (knots) , 

W500 = 500-mb wind speed (knots) . 

The SWEAT index is always positive. No individual term may 
ever be negative. Based on empirical data, the SWEAT index threshold 
value for tornado cases is ~350, while for severe thunderstorms it 
is ~250. Miller [35] refers to this initial SWEAT index as the 
"Soft SWEAT" index. 

The SWEAT index was further modified [32] to include the 
500-mb/850-mb level wind directional shear term. This shear term is 
also based upon directional wind shears observed during severe 
weather cases and changes the SWEAT index (SWEAT^) equation to read: 

SWEAT 2 = 12 T d85Q + 20 (TTI - 49) +2W 85Q + W 5Q() + 125(S+0.2) , (2) 


I I II III I 
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where 


S - sin(WD 5Q0 - WD g50 ) , 
and 

WD = wind direction (degrees) . 

If the 850-mb wind is not within the range 130 and 250 deg, 
or if the 500-mb wind is not between 210 and 310 deg, or if the 
expression WD 50Q - WD g50 < 0, set. S = -0.2 to drop the shear term. 

The addition of the shear term to the SWEAT index raises the severe 
thunderstorm threshold to ~300, and that for tornadoes to ~400 to 425. 

The SWEAT index is not a tool for forecasting ordinary 
thunderstorms. It is designed to indicate the potential of severe 
thunderstorms (with gusts at least 50 kts and/or hail at least 
0.75 in. diameter) or tornadoes. 

Lastly, in the SWEAT equation. Miller [35] replaced the 
850-mb level with the 900-meter level (except in the TTI and shear 
calculations) and changed the wind directional shear procedure. The 
revised SWEAT equation (SWEAT g ) thus reads: 

SWEAT 3 = 12D e + 20(TTI - 49) + 2W g + W 500 + 125 f(a) , (3) 
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where. 


n = low level dew point (°C) at the 900-meter level J 
e 

2 

= low level wind speed (kts) at the 900-meter level, 

3 4 

f (a) = a step function of the veering angle W g to W 5Q0 . 

The term f(a) is set to 0 if both W e and W 5Q0 are not >15 kts. The 
850- and 500-mb level wind directions must also fall within earlier 
stated ranges (see Figure 13). All other terms are defined exactly 
as before. The use of 900-meter level parameters in the calculation 
of the SWEAT index is referred to as "BLM SWEAT" since it is the 
equation used in the AFGWC Fine Mesh and Boundary Layer Models (BLM) 
forecast model . 

The soft SWEAT index plots can be computer-calculated within 
1.5 hours of the 00 GMT or 12 GMT sounding time. The BLM SWEAT 
calculations take up to four hours of computer time. Both SWEAT 
index methods are currently being used and 1 2- , 24-, and 36-ho'ur 
SWEAT index prognostic maps are generally output. 

Recently, Miller and his associates [35, 36] have noticed 
that many times severe weather has formed within overlapping areas 


^Use 850-mb dew point in soft SWEAT. 

^Use 850-mb wind speed in soft SWEAT. 

^Use of the sine function was discontinued for soft or BLM 
SWEAT because it was not representative from 30 to 120 deg. 

^Veering is defined as a change in wind direction versus 
altitude, in a clockwise sense. 
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0 <20 <40 <60 <80 <100 <120 <140 <160 <180 

a = Veering Angle (degrees) 

Figure 13. Step function used in computation of SWEAT 
veering term. 
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of high SWEAT and high SPOT (Surface Potential index [37]) values. 
Therefore, these two statistically derived indices can be used 
together as an aid to accurately identify short-term (three to six 
hours), small-scale potential severe storm areas. The false alarm 
rate within the SWEAT/SPOT forecast overlap area is much smaller 
than that of either index used separately. 

The second SWEAT equation (Eq. (2)) has been programmed and 
is used in the present study as the SWEAT index. 

Vertical -T otal s I ndex 

In 1967, Miller [38] introduced the term "vertical totals" 
in relation to potential thunderstorm development. The vertical- 
totals index (VTI) represents the stability of the atmosphere 
(temperature lapse rate) between 850 and 500 mb with no moisture 
parameters involved. It is defined as the 500-mb temperature sub- 
tracted from the 850-mb temperature; that is, 

VTI = t 850 " T 500 ^ 

Vertical -total s values give a measure of instability. 
Generally, values >26 represent thunderstorm development without 
regard to moisture. Specific areas and their approximate critical 
VTI thunderstorm threshold values are listed in Table 13. 

Although the VTI can be used alone, it is also valuable when 
added to the cross-totals moisture index (CTI). Combination of VTI 
and CTI results in a total-totals index (TTI) is described later. 
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Table 13. Vertical-Totals Index Thunderstorm 
Threshold Values for Different Areas 


Area 

Critical VTI 

Gulf Coast 

>26 

British Isles 

>22 

Western Europe 

>28 

West of the Rockies 

>29 

Pacific Coastal Areas 

>30 

Great Lakes 

>30 
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Cross-Totals Index 


Also in 1967, Miller [38] introduced the cross-totals index 
(CTI) as the 500-mb temperature subtracted from the 850-mb dew-point 
temperature; that is, 

CTI = T d850 ' T 500 ^ 

Thus, a low-level moisture parameter is introduced into the index 
calculation. The CTI has been used to indicate thunderstorm 
potential, with the cross-totals thunderstorm threshold usually 
about 18. However, along the Gulf Coast a CTI of 16 (with VTI >23) 
generally produces a thunderstorm. The cross-totals index is also 
an initial index used in the calculation of the total-totals 
stability index explained in the next section. 

Total-Tot a ls Index 

In 1967, Miller [38] introduced the concept of the total - 
totals index (TTI) as being a measure of atmospheric instability 
between the 850- and 500-mb level. The TTI is defined as the 
arithmetic sum of the vertical -total s index and the cross-totals 
index; that is, 


TTI 

= VTI + CTI (°C) , 

(6) 

TTI 

(T 850 + T d850^ ' 2 T 500 ( ° C) * 

(7) 


The VTI thunderstorm threshold of 26 and the CTI of 18 pro- 
duces a minimum threshold of 44 for the total -totals index. Total - 
totals index values >50 generally indicate the potential of numerous 
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and severe thunderstorm/tornadic activity if an adequate low-level 
moisture supply and a trigger mechanism are both present. The TTI 
proved to be more accurate in forecasting of thunderstorms, in all 
places and seasons, than did either the VTI or CTI alone. 

Theta-E Methods 

The use of equivalent potential temperature (0 £ ) can be used 
in synoptic meteorological practice as a measure of atmospheric 
stability [39]. The quantity 0^ is quasi -invari ant (conservative) 
with respect to both dry and moist adiabatic processes, and is 
invariant (does not change) with respect to evaporation of falling 
rain [40, 41, 42]. Equivalent potential temperature is a single 
parameter which takes into account both temperature and moisture 
content. Theta-E cannot be measured directly since it is simply a 
concept. It is defined as follows: A parcel of air at tempera- 

ture T , dew-point temperature T^ , and any pressure level P q rises 
vertically by a dry-adiabatic process until saturated (at LCL) and 
then follows the moist-adiabat until all moisture precipitates out. 

At this point, the moist-adiabat is parallel with the dry-adiabat on 
the Skew-T diagram. If the air parcel is now compressed dry- 
adiabatical ly down to a pressure of 1,000 mb, it will have a tempera- 
ture defined as the equivalent potential temperature expressed in 
degrees absolute. Figure 14 illustrates this process. 

Theta-E is also a measure of potential stability in that it 
gives a measure of the effect lifting will have on a column of 
air [39]. Theta-E can be computed at two vertical levels on a 
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sounding, and if it decreases with height between the two levels 
(i.e., AG^/AZ or AQ^/AP < 0), this layer is absolutely unstable if 
lifted to the saturation level [26, 31]. On the other hand, the 
layer remains stable if lifted to the saturation level when 0^ 
increases with height. Reference [43] reports that computed 700-mb 
Theta-E charts were being transmitted via the facsimile network to 
aid in the forecasting of thunderstorm activity as early as 1950. 

The 700-mb Theta-E critical value of ~327°K together with the 6 g/kg 
mixing ratio line was generally used to outline areas likely to 
experience heat-type thunderstorms. Values of 0^ ~321°K and 
w~4 g/kg indicated the potential of a lifting-type thunderstorm. 

Delta Theta-E (A6 ) values have also been used in thunder- 
storm forecasting [44], which expressed the change in 6^ versus 
pressure-altitude (AG^/AP), as indicated earlier. The difference in 
Theta-E between 850 and 700 mb forms a lower index, and that between 
700 and 500 mb an upper index; that is. 


II 

1 

LU 

CD 

< 

e E850 " e E700 ’ 


(8) 

A0 EU - 

e E700 " e E500 ' 


(9) 

i index 

are presented in Table 14, 

with 


licating 

instabi 1 ity. 



et al . 

[31] have used and tested 

a simple 

e E 


difference equation of the form: 


0 * = 0 


e ESfc + e E850 


E700 


(10) 
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Table 14. Delta Theta-E Upper and Lower Index Critical Values 


A9 E Index Stabi 1 ity 


> ±0 

Lower 1 

Unstable 

> -5 

Upper J 


< ±0 but > -2 

Lower 1 

Questionably 

Unstable 

> -5 

Upper J 


< -2 

Lower 1 

Stable 

< -5 

Upper | 
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where 0* defines convective instability, if 0* < 0 at altitudes 
between 700 mb and close to the ground (surface and 850 mb). This 
difference is similar to the lower index of the Delta Theta-E method 
mentioned previously. The 0^ index selected for this study is that 
of Alaka, with Eq. (10) being modified by replacing 0^ c with 
since the 900-mb level is the lowest level of averaged data available 
This index will be referred to as 6*. Also, all equivalent 
potential temperature (0^) computations made in this study are 
derived from the approximate form (0^) from Eq. (21), as explained 
later on in the Energy Index section of this chapter. 

Showalter Stability Index 

The Showalter Stability Index (SI) was developed by 
A. K. Showalter of the United States Weather Bureau in 1946 and 
documented more widely in 1953 [5]. It is a thermodynamic static 
index which can provide a quick, simple estimate of possible 
thunderstorms based on the potential (convective) instability concept 
This index was designed for initial use in the southwestern United 
States, but has been used extensively all around the world. It was 
believed that areas of instability are not generally altered signifi- 
cantly at 850 mb and above. Therefore, a stability-index map based 
on this level and above can be derived and the stability movement 
prognosticated for 12 and up to 24 hours. 

The SI is computed as follows: Dew-point temperature (T^) 

and temperature (T) values (in °C) are obtained at the 850-mb level 
(the assumed top of the moisture layer), together with the tempera- 
ture value (in °C) at 500 mb. Showalter indicated that mountain 
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sites can use T and T. values from the 700-mb level instead of the 

d 

850-mb level in their SI computations. The 850-mb parcel is now 
lifted dry-adiabatically to the saturation level (LCL) and then 
lifted moist-adiabatical ly to 500 mb. The lifted 500-mb temperature 
is then subtracted from the observed 500-mb temperature; that is,. 


SI 


t obs 

500 


VlFTED 

500 


(in °C) . 


cm 


The procedure is depicted graphically in Figure 15. Showalter 
Stability Index values of +3 deg or less generally indicate probable 
showers and some thunderstorms in the area; SI values from +1 to 
-2 deg indicate increasing probability of thunderstorms; SI values 
from -3 to -5 deg (or less) indicate possible severe thunderstorms; 
SI values from -6 deg or less indicate suspect conditions for 
tornadoes. 

The Showalter index has been used extensively over the years 
in different capacities. It has been directly correlated with 
hail [45, 46] and with storm radar echoes [47 through 49]. It has 
been used in the forecasting of general showers resulting from 
surface heating as well as from lifting [3]. This index has 
also been used in heating calculations because, besides being a 
function of the 6^ or 0^ lapse rates, it is also partly a 
function of the ordinary temperature lapse rate and is, therefore, 
indicative of stability for use in surface-parcel heating appli- 
cations. This index is a measure of convective stability when the 
index value is greater than +6, and convective instability when 
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Figure 15. Showalter stability index computation method. 
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values are less than zero. Also, the index is a first approximation 
in estimating latent instability, because negative index values do 
indicate that a positive area (energy) does exist above the LFC. 

The Showalter index is limited for use in mountain areas, 
and will not work well if the air is extremely dry, or if critical 
instability exists higher than the 850-mb level. This is because 
the index uses only the one lower point at 850 mb as being repre- 
sentative of low-level moisture and temperature. 

Rackcliff Instability Index 

In 1962, Rackcliff [50] introduced a simple latent insta- 
bility index, patterned after the lifted-index [51], for use in 
regional forecasting of air-mass-type summer thunderstorms in the 
British Isles and Western Europe. 

While the lifted-index uses a forecasted maximum afternoon 
temperature in its calculation, Rackcliff used a computed tempera- 
ture in the calculation of his index. The 900-mb wet-bulb potential 
temperature ( 6 W 9 qq) was the low-level temperature parameter selected 
by Rackcliff. It is obtained by taking the 900-mb wet-bulb 
temperature and descending moist-adiabatical ly to the 1,000-mb level, 
as shown in Figure 16. The e w gQQ value is believed to be repre- 
sentative of the air at low levels and is also only slightly 
affected at night by outgoing terrestrial radiation. The environ- 
mental temperature at 500 mb (T^qq) is a 9 a in used as the indicator 
of middle-tropospheric temperature. The Rackcliff index (RI) is 
then defined as the algebraic difference of the 500-mb temperature 
from the 900-mb wet-bulb potential temperature; that is. 
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Figure 16. Rackcliff instability index computation method. 


58 


where positive values represent latent instability. Rackcliff 
determined the following thunderstorm/no thunderstorm criteria: 


1. RI < 25 (stable condition). 

2. RI > 25 (showers possible). 

3. RI > 30 (thunderstorms possible). 

4. RI > 35 (heavy thunderstorms possible). 

The value of 30 is a thunderstorm threshold value used in fore- 
casting nonfrontal thunderstorm activity in the British Isles. 

Jefferson Instability Index 

A modification of Rackcliff's index was made by Jefferson [52] 
in 1963 so that the instability index could be used in summertime 
air-mass thunderstorm forecasting at the London Airport. Jefferson 
determined that Rackcliff's index makes no allowance for the fact 
that instability in a layer depends not only on the temperature 
difference across the layer, but also on its mean temperature. Since 
the value of 0^ varies between 10 and 20°C over northwest Europe in 
summertime thunderstorm situations, this would give a variable 
Rackcliff index value between 36 and 29. Therefore, Jefferson 
amended Rackcliff's formula with an empirical study and obtained an 
instability index value independent of temperature, but with the 
same threshold value of 30 for thunderstorms. This was true for a 
wide range of temperatures. This modified index can now be used in 
wider areas and for all seasons. The Jefferson instability 
index (JI) is expressed as: 


JI ' 1 - 6 e W 9oo " T 500 11 » 


03 ) 


where ^goo is the wet "bulb potential temperature (°C) and 

T 500 is the observed 500-mb temperature (°C). Positive index values 
represent instability. 

Modified Jefferson Instability Index 

In 1963, Jefferson [53] published a second modification to 
the Rackcliff index, or simply, a modified Jefferson index (MJI). 
While using the Jefferson index at the London Airport, it was 
determined empirically that the index was forecasting thunderstorms 
(i.e., JI values exceeded 30) in the Mediterranean area, but many 
times no thunderstorms formed. This was found to be caused by very 
dry air existing above 900 and 500 mb over the Mediterranean area. 
Since the base of thunderclouds over the Mediterranean is generally 
quite high (-700 mb), the idea of introducing a 700-mb moisture 
parameter seemed logical, as long as the index continued to work for 
north-central Europe. This modified Jefferson index (MJI) is 
written as: 

MJI = 1.6 e w goo ' T 500 " 1/2 AT d700 ‘ 8 ( c ) » ( 14 ) 

where, 

0 W9OO = wet-bulb potential temperature (°C) , 

T cnn = 500-mb observed temperature (°C) , 
bUU 

AT d700 = 7 00-mb dew-point depression (T 70Q - T d70Q ) (°C) . 
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The factor 1/2 ATjygg was introduced to avoid overweighting by the 
AT d7 go parameter. 

Boyden Instability I ndex 

Just prior to the 1963 publication of the modified Jefferson 
index, Boyden [54] also introduced an instability index to be used 
in the forecasting of thunderstorms and heavy rain over southeast 
England during the months of May to September. Boyden assumed that 
the development of heavy showers and thunderstorms over land on a 
summer afternoon depends on the mean temperature lapse rate only up 
to 700 mb. For neutral static stability conditions (i.e., dry-bulb 
temperatures along a moist adiabat), Boyden determined that the 
1,000- to 700-mb thickness (in decameters) minus the 700-mb 
temperature (°C) was an approximate constant (~294) for all summer- 
time atmospheric conditions measured over Crawley, England. 
Instability is then measured by the amount this difference exceeds 
the constant. Therefore, instability exists if the 700-mb tempera- 
ture is a low (cold) value as compared to the 1,000- to 700-mb 
thickness value. Boyden 's index (BI) is expressed as: 

BI = AZ ( 1,000 to 700) ‘ T 700 " 200 * < 15 

where, 

AZ = 1,000- to 700-mb thickness (decameters) , 

T 700 = 700-mb temperature (°c) . 

The units conflict in the BI expression. Only the numerical value 
should be used. The value 200 is used to remove the large unwanted 
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number generated by this index. It allows the BI to take on a value 
around 90. 

The BI is strictly a measure of the mean stability in the 
entire layer below 700 mb. The Boyden index is not intended to 

r 

forecast slight or moderate showers. The diurnal variation of BI 
was found to be low, allowing a 12-hour forecasting of the index to 
be made. It was determined that Boyden index isopleths (drawn in 
intervals of two units) move with the 700-mb wind. 

For both frontal and non-frontal summer days, it was found 
that there was, indeed, a marked increase in thunderstorm/heavy rain 
occurrence when BI reached values of 94 and higher. Since humidity 
was found to be very loosely related to the development of thunder- 
storms, it was not included with the Boyden index. The main 
advantage claimed for the Boyden index is its usefulness at mobile 
sites during frontal or non-frontal weather. 

For the present study, the Boyden index was modified, since 
1,000-mb heights are not obtainable from the averaged soundings. 
Therefore, the 900-mb height was used in place of the 1,000-mb height 

K-Index 

The K-index (KI) was developed by Whiting and documented by 
George (both of Eastern Air Lines) in 1960 [55]. This simply derived 
stability index is used in the forecasting of inland air mass thunder 
storms with weak winds and without apparent frontal or cyclonic 
influence. It is prepared from the 1200 GMT soundings and is 
generally issued on an areal map (with KI intervals every five units) 
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The Whiting-George K-index measures air mass thunderstorm potential 
by direct indication of the vertical temperature lapse rate 
(Tggo - T 5 qo^’ ^ ower atmospheric moisture (T^SSO^’ and very indirect 
indications of the vertical extent of the moist layer (700-mb dew- 
point spread). The K-index is expressed as: 

KI = ( T 850 " T 500^ + T d850 " ^ T 700 " T d700^ ^ 

where K-values versus thunderstorm occurrence frequencies generally 
fall within the categories given in Table 15. 

The K-index map used concurrently with a subjective analysis 
of convergence and relative vorticity has been proven by George to 
be a valuable air mass thunderstorm forecasting tool. Areas of 
confluence, determined by constructing 850- plus 700-mb height areal 
charts, are used to represent convergent flow conditions between 
these two levels. Confluence areas below 700 mb, with winds 
<20 knots, generally require an adjustment to the next higher 
category of K-values. If the winds are >20 knots, adjust upward two 
categories. Positive vorticity also increases the chance of thunder- 
storm development. 

Bryan [56] and Hambridge [57] have tested the K-index versus 
thunderstorm activity over the mid-South and Western United States, 
respectively; they found a high correlation. Hambridge suggested 
the assignment of thunderstorm probabilities versus K-value given in 
Table 16. 

In 1971 [58], the K-index chart was added to the lifted 
index panel of the composite moisture index chart. This chart is 
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Table 15. K-Index Thunderstorm Threshold Values 


K-Index Value 

Thunderstorm Frequency 

K < 20 

None 

20 < K < 25 

Isolated 

25 < K < 30 

Widely Scattered 

30 < K < 35 

Scattered 

35 < K 

Numerous 
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Table 16. K-Index Thunderstorm Probabilities 


K-Index Value Thunderstorm Probability 



< 

15 


0% 

15 

to 

20 


< 

20% 

21 

to 

25 

20 

to 

40% 

26 

to 

30 

40 

to 

60% 

31 

to 

35 

60 

to 

o 

00 

36 

to 

40 

80 

to 

90% 


> 

40 

Near 

100% 
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distributed via the NWS NAFAX (National Weather Service National 
Facsimile Network) system to all meteorologists across the United 
States. 

Bradbury Potential Stability Index 

In 1977, Bradbury published an article [59] dealing with the 
use of wet-bulb potential temperature (0^) charts in weather 
analysis and forecasting. One conclusion he reached was that many 
summer thunderstorms broke out over Europe when low-level southerly 
winds advected air with ^ggg > 16°C. Bradbury then developed a 
potential stability index (BPI), since the Q^q values alone failed 
to identify occasions of thunderstorm development in a relatively 
cool air mass. This index, similar in structure to Rackcliff's and 
Jefferson's index, is defined by subtracting the value of 6^ at 
850 mb from the value at 500 mb. A negative value of this difference 
indicates that the air between the two levels is potentially unstable. 
In equation form, the BPI is expressed as: 

BPI = 0 W5OO " 0 W85O ^ C ) * ( 17 ) 

One can obtain the BPI from a thermodynamic diagram procedure as 
illustrated in Figure 17. 

Bradbury also found that the BPI varied as a function of 
0 W85O’ w ^ en usec * 1n t * ie f° recas t'’ n 9 of thunderstorms during the 
year. This is illustrated in Figure 18, where 5%, 50%, and Limit 
represent the cumulative percentage frequency of BPI versus e ^ggQ 
for 544 thunderstorm day soundings from 1973 to 1976. The graph 
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Figure 17. Bradbury potential stability index computation method. 



0 2 A 6 8 10 12 M 14 18 20 °C 


e W850 mb 

Figure 18. Relationship between Guocq and the BPI 
on thunderstorm days [59]. 
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merely indicates a range of conditions that existed when thunder- 
storms occurred. Thunderstorms would be unlikely outside the range 
given. One should not use the BPI as a strict thunderstorm fore- 
casting rule, but rather, as a guide along with routine surface and 
upper-air charts. 

Energy Index 

A unique and reliable substitute for the widely used thermo- 
dynamic indices, used in the forecasting of convective storms, is 
the total energy index (El). It was introduced by Darkow [60] in 
1967 and deals with the total energy (Ey) of a unit mass of air. 

The specific enthalpy (c T), potential energy (gZ), latent 

2 

energy (Lq), and kinetic energy (W /2) of the unit mass of air is 
combined as: 

E T = c p T + gZ + Lq + W 2 /2 (cal gm -1 ) , (18) 

where, 

Cp = specific heat of air (0.24 cal gm ^ °K ^ ) , 

T = temperature (°K) , 

_2 

g = acceleration of gravity (980 cm sec ) , 

Z = altitude (km) , 

L = latent heat (cal gm ^ ) , 
q = specific humidity (gm kg~^) , 

W = scalar velocity (cm sec ^) . 
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Since the kinetic energy term is two orders of magnitude smaller 
than the other three terms, it can be neglected, resulting in the 
energy formula being called static energy (E s ); that is, 

E s s Ej » CpT + gZ + L q w , (19) 

where, 

L ~ L Q = 600 cal gm -1 , 
q z w = mixing ratio (gm kg' 1 ) ; 

therefore, 

E t s 0.24 T + 2.34 Z + 0.6 w (cal gm' 1 ) . (20) 

Static or total energy is conserved with respect to both 
types of adiabatic processes and is related to the pseudo-equivalent 
potential temperature (0^) and wet-bulb potential temperature (0 ). 
This fact can be seen by dividing E T by c^, which produces a geo- 
equivalent potential temperature (0,,.-), which is a conservative 

bt 

(invariant) property in regard to adiabatic processes. The term 9g^ 
differs just slightly in definition from 0^ and is expressed as: 

E t 

6 = — = T + 9.8 Z + 2.5 w (°K) . (21) 

GE c p 

Total energy or geo-equivalent potential temperature can 
both be computed easily for use in the forecasting of convective 
activity. This total energy concept can be used in both ascent and 
descent air parcel theory convective calculations, and the amount of 
potential convective instability of the air column is indicated by 


69 



the decrease of total energy with increasing altitude. This defines 
the Darkow total energy index (El). It is expressed as the algebraic 
difference between the atmospheric total energy at the 500- and 
850-mb levels; that is, 

EI = E T500 _ E T850 (cal gm' 1 ) . (22 ^ 

Empirical testing of the index produced the ranges given in Table 17 
for forecasting severe weather. The total energy index horizontal 
map pattern turns out to be very similar in structure to the 
Showalter index pattern. This is due in part because low-level 
total energy is usually greater than mid-tropospheric values. 

The total energy index combines temperature, moisture, and 
height fields. Darkow [60] indicates that this gives it a possible 
advantage over the Showalter and lifted indices since it is the only 
one to take into account the possible contribution of descending, 
potentially cold, mid-tropospheric air on the total energy release 
of convective storms. Most indices involve only the process of 
ascending warm air. 

Darkow took an additional step by suggesting that a modified 
energy index can be developed which takes into account the mean 
mixing ratio of the lowest 100-mb layer (or of the first kilometer 
altitude) above the ground. This may be more representative of 
lower level moisture than using just the 850-mb value of mixing 
ratio. 

A number of atmospheric studies have used the Darkow energy 
index and total energy concept. Some of these studies are presented 
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Table 17. 


Energy Index Values Used in Convective 
Forecasting 



El 

Forecast 3 

0.0 

to -1.0 

Non-severe thunderstorms 
possible. 

-1.0 

to -2.0 

Isolated severe thunder- 
storms possible. 


< -2.0 

Severe thunderstorms and 
tornado activity possible. 


a If a trigger mechanism is available to 
release potential instability, otherwise, convective 
activity may not take place. 


in References [61 through 64]. Eagleman [61] used Darkow's index 
separately and also combined with a wind shear index to aid in fore- 
casting tornadoes. Darkow has also applied the static energy concept 
to surface analysis in detecting areas of high static energy related 
to thunderstorm and severe storm occurrences [65, 66]. 

Martin Index 

Wright-Patterson Air Force Base, Ohio, published a stability 
index constructed by D. 0. Martin [67] said to be more sensitive to 
low-level moisture than the Showalter index, since it uses the 
maximum value of low-level moisture. 

The procedure in computing the Martin index (MI) (see 
Figure 19) is as follows: From the 500-mb temprature (A), descend 

moist-adiabatical ly to the intersection (B) of this line and the 
mixing ratio line that passes through the point (C) of maximum 
mixing ratio. From this intersection, move dry-adiabatically to the 
850-mb level (D). The MI is defined as the difference between the 
observed sounding temperature and calculated temperature at 850 mb; 
that is, 

MI = T 850 - T 850 ( ° C) • (23) 

Calc. Obs. 

The only exception to this procedure occurs whenever a marked 
low-level turbulence or subsidence inversion (non-surface, non- 
radiation) is established below 850 mb. Then the point (D) is 
obtained at the pressure level where the inversion base is located. 

The normal and the exception cases are illustrated in Figure 19. 
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Figure 19. Martin stability index computation methods. 


Since this study involves the 900-mb level as being the 
closest level to the surface, the index will be referred to as the 
modified Martin index (MI). 
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CHAPTER V 


AVE-IV PROFILE ANALYSES 
A. Introduction 

Before stability indices can be constructed, used, or evalu- 
ated, the atmospheric parametric profiles themselves need to be 
examined and understood. Therefore, this section presents a dis- 
cussion of the AVE-IV average and lag soundings as they are compared 
with the four MDR precipitation categories, and as they are compared 
with each other. Tabular values of these profile parameters have 
been presented in Chapter III, Tables 4 through 11, pages 32 
through 39. 

For clarity, and to avoid confusion, the average profiles 
which pertain to precipitation conditions occurring at the time of 
the sounding observation will hereafter be referred to as AVG. Also, 
the averaged lag profiles, which represent the environmental obser- 
vations three hours prior to a precipitation category occurrence, 
will hereafter be referred to as LAG profiles. 

Throughout this section, more attention will be given to the 
parametric average profile differences which exist between precipi- 
tation categories A and D (defined in Table 3, page 27). 

Categories B and C parameter differences have been compared but are 
not always presented here because they either do not represent any 
drastic environmental change, or since they do have inherent 
category D information, they could present a bias. In most all 
cases these two intermediate precipitation categories merely link 
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categories A and D. The average profiles of LAG and AVG temperature, 
potential temperature, and winds are presented and compared for the 
storm categories A and D. Differences are noted for possible 
inclusion in a forecast-type of storm index. 

B. AVG Profile Comparison 


Temperature/Moisture 

The vertical temperature profile differences noted between A 
(non-precipitation) and D (severe storm) conditions for the AVE-IV 
average (AVG) profiles are shown in Figure 20. Although the two 
temperature profiles are almost identical above 700 mb, the D profile 
temperatures increase 5°C warmer than the A conditions between the 
700- and 900-mb levels. 

The dew-point temperature profile for D conditions is 
2 to 10 deg warmer (more moist) than for A conditions at all alti- 
tudes, as shown in Figure 20. Most of the difference (6 to 10°C) 
occurs between the 600- and 800-mb levels. 

Winds 

As one would expect, for all altitudes, winds are higher 
when going from precipitation categories A to D using the AVG wind 
profile information. Wilson and Scoggins [20] also confirmed this. 
This increase pertains to both meridional and zonal wind components. 
The category D meridional wind component exhibited the most 
difference (~8 m sec" 1 ) over category A conditions. Zonal (and 
scalar) wind differences between A and D categories generally range 
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Figure 20. AVG temperature and dew-point profiles for A and D MDR conditions. 




between 2 and 6 m sec - ^ (see Figure 21). All wind components calcu- 
lated are positive ( i . e . , zonal winds being westerly and meridional 
winds southerly) for all AVG and LAG conditions presented. Zonal 
westerly winds dominate in magnitude. 

C. LAG Profile Comparison 

To determine if a forecast scheme can be realized based on 
AVE-IV data, this stability study will involve the analysis of LAG 
profiles and how their average conditions differ from the AVG 
profiles representing storm activity. 

Temperature/Moisture 

LAG thermodynamic profile conditions for the four precipi- 
tation categories are very similar in appearance to the four 
respective AVG profiles. Temperatures of category D are warmer by a 
similar magnitude than category A, as was the case for the AVG pro- 
files. This effect extends higher, however, from 900 to 650 mb. 

The category D dew-point profile also remains warmer by a similar 
moisture differential spread, as was the case with the AVG dew-point 
data. 

Winds 

LAG winds again invoked a pattern similar to that of AVG 
winds, with both D wind components exhibiting stronger flow than 
A wind component conditions. Category D LAG wind directions between 
7 and 13 km altitude are slightly westerly, so as to resemble 
category B LAG wind directions over this altitude range. This is 
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Figure 21. AVG wind components for MDR A and D conditions. 




completely unlike the AVG D, AVE-IV storm wind directions, which are 
farthest away from the west of all four precipitation categories. 
This directional change will be discussed further in Section D. 

D. AVG/LAG Profile Comparison 


Winds 

Differences between the AVE-IV average LAG (three to six 
hours prior to storm occurrence) and AVG (time of storm occurrence) 
wind profiles are, again, generally small. However, the largest 
differences do occur between the category D (severe storm) wind 
profiles of each. Therefore, only category D comparisons will be 
discussed here. 

From Tables 7 and 11, pages 35 and 39, respectively, 
scalar wind speeds and U-component (zonal) speeds are ~2 m sec"^ 
stronger for the LAG average than for the AVG average, between 
6- and 12-km altitude. Wind magnitude differences were less than 
this value above 12-km altitude (see Figure 22). 

Magnitudes of the V-component wind (meridional) give slightly 
stronger (~2 m sec - ^) southerly winds at AVG time than LAG time, 
between 8- and 12-km altitude. This stronger AVG V-component effect 
coupled with the weaker AVG U-component results in the AVG wind 
direction between 8- and 12-km altitude being ~5 deg more from the 
south (248 deg) than the LAG average (253 deg). Meaning, on the 
average, winds during the LAG period are 1 to 2 m sec”^ stronger and 
from a more westerly direction than conditions existing during severe 
storm occurrence. 
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Figure 22. Wind components for AVG and LAG type D MDR conditions. 


Temperature/Moi sture 


The thermodynamic structure between AVG and LAG category D 
profiles, in terms of temperature and dew point, is shown in 
Figure 23. Average conditions for each respective parameter are, 
indeed, very similar. The unusual feature is the ~1.5°C temperature 
difference that exists around 650 mb, with equal or lesser differ- 
ences indicated between the 800- and 500-mb levels, and with LAG 
temperatures being slightly warmer. Also, the layer between 500 and 
650 mb is more unstable three hours prior to storm activity. The 
temperature lapse rate of the LAG sounding between these two levels 
is 1 5 . 4°C/1 50 mb; whereas, only 14.4°C/150 mb (difference = 1.0°C) 
existed during AVG storm time. This slightly more unstable layer is 
noticed at a lower altitude between 800 and 650 mb at storm time. 

It then has a temperature gradient of l 2 . 0°C/ 1 50 mb as compared to 
1 0 . 9°C/1 50 mb (difference = 1.1°C) for this layer on the prior LAG 
profile. 

Three hours prior to storm activity, the dew-point temperature 
at 900 mb is ~1°C higher than at the time of storm activity. By the 
850-mb level, this difference vanishes and neither AVG nor LAG dew- 
point temperatures dominate above this level. 

E. Theta-E AVG/LAG Comparison 

Before the stability index results are presented and dis- 
cussed, it is desirable to select an index or procedure involving 
equivalent potential temperature (0 £ ) as an instabi 1 ity measure (see 


82 



PRESSURE (mb) 







Chapter IV discussion of 0 £ ). This is an important decision to make 
because of the many number of ways in which 0^ can represent atmos- 
pheric instability. 

AVG 0 ^ profiles are presented for the four MDR conditions in 
Figure 24. Similar LAG 0^ profiles are given in Figure 25. The two 
figures present similar results and show higher values of 0^ for more 
severe weather activity. The comparison between LAG and AVG profiles 
of 0£ are given in Figure 26 for categories A and D MDR conditions. 
The LAG profiles exhibit slightly greater 0^ (or total energy) than 
do the AVG profiles. The altitude of minimum 0^ occurs at ~700 mb 
for the two profiles with no MDR activity (category A), while it 
occurs higher ("“600 mb) for both the MDR>7 (category D) profiles. 

One item of significance is the more stable 0^ gradient observed 
between 850 and 800 mb on the AVG-D profile than the LAG profile 
indicates. The AVG-D 0^ profile also indicates a slightly more 
unstable region between 750 and 700 mb, as compared to LAG-D 
conditions . 

Results from the convective stability equation of 
Alaka et al . [31] and the 0^ differences between 800 to 850-, 

800 to 900-, and 700 to 750-mb levels, as suggested by observing the 
0 E LAG and AVG vertical profiles, are presented in Table 18. It 
should be noted that Alaka's equation (Eq. (10)) is similar in 
structure to the Delta Theta-E equation (Eq. (8)), as discussed 
earl ier in Chapter IV. 

While the 0* equation of Alaka (Eq. (10)) is representative 
of the entire lower atmospheric instability when applied to 
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Figure 26, LAG and AVG equivalent potential temperature (0 F ) vertical 
profiles for type A and D MDR categories. L 


Table 18. Theta-E Differences (° 


LAG 


Pressure Level MDR Conditions 

Difference Category ABC 


00 

00 

1. 

Alaka Bjr (700 to surf.) 

-3.47 

-3.17 

-6.02 


2. 

A0 e (800 to 850) 

-1.15 

-0.85 

-1.05 


3. 

A6 e (800 to 900) 

-2.25 

-1.34 

-2.14 


4. 

A0 e (700 to 750) 

-0.36 

-1.17 

-2.57 



individual soundings, it is not when applied to averaged soundings, 
as shown in Table 18. The layer between the surface and 700 mb 
appears to be too large, to note meaningful differences (only 7%) 
between LAG-D and AVG-D conditions. 

Values of A0 £ over narrower layers are presented in Table 18, 
using two lower atmospheric levels (items 2 and 3 of Table 18) and 
one upper atmospheric level (item 4). The significant point to 
notice is that pre-storm conditions (LAG-D) have a stronger 
0£ lapse rate between 800 mb and below, as compared to storm con- 
ditions (AVG-D); the AVG being 38 and 53% lower than LAG conditions. 
The opposite is true when the 700 to 750-mb 0^ gradients, for both 
D MDR conditions, are compared. Here, AVG-D gradients are 49% higher 
than LAG-D conditions. 

The results presented in Table 18 do indicate that possibly 
a stability index could be arrived at by using one or two A0^ 
parameters in its computation. This would make the index a function 
of both a temperature and moisture input over a 50- to 100-mb spacing 
in the atmosphere. 

F. AVG/LAG Conclusions 

A few general conclusions can be made regarding the atmos- 
pheric environment three hours prior to severe storm development and 
during the occurrence of the severe storm. They are as follows: 

1. The LAG profile exhibits a more unstable temperature 

gradient (15. 4°C/1 50 mb) in the upper atmosphere between 
650 and 500 mb three hours before severe storm activity. 
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By the time of severe storm occurrence, this unstable 
layer has fallen 150 mb and is located between 650 and 
800 mb with a 12.0°C/150 mb gradient. 

2. Dew-point temperatures are warmer (more humid) from 900 
to 850 mb by 1°C three hours prior to severe weather. 

3. Scalar wind and zonal wind speeds are stronger by 
~2 m sec”^ three hours prior to storms. At storm 
occurrence time the 2 m sec -1 stronger, southerly, AVG 
meridional wind component, located between 8 and 12 km 
altitude (400 to 200 mb), results in producing a more 
southerly direction (248 deg) than three hours prior 
(253 deg). 

4. Equivalent potential temperature (0^) differences 
between 800 and 850 mb are the most unstable three hours 
prior to convective weather; whereas, 0^ differences 
between 700 and 750 mb are more unstable at the time of 
severe weather occurrence. 

The results obtained here could be used in the construction 
of a type of three-hour lead time severe storm index. However, 
these changes noticed in the atmospheric structure are, indeed, all 
very small changes. It should be noted that they are small due to 
the fact that they are based on the average of many pre-storm 
soundings taken during only one independent major storm system 
development/movement. It may also turn out that these AVE-IV con- 
clusions may or may not apply to a different storm sounding history 



for a station. One purpose of the present study is to test this 
theory, thereby affecting current understanding of the severe storm 
environment. 


G. Possible Stability Index 

At this point, a stability-type of index could be constructed 
consisting of a temperature difference (AT), or an equivalent 
potential temperature difference (A0^), between two levels. A low- 
level dew-point indicator (AT^), and possibly a wind magnitude term 
(AW) and directional term (AWD) could be included. The combination 
of these terms, with the appropriate multiplication weighting 
factors (M), could result in a meaningful severe storms lag index 
(SSL). A possible form of the equation is: 

SSL = M-| (AT ) + M 2 (AT d ) + M 3 (AW) + M 4 (AWD) + M 5 (A0 E ) . (24) 

All or only a couple of the terms expressed in Eq. (24) may prove 
useful as an index parameter when compared to its respective severe 
storm threshold value. Considering the broad scale in which these 
five terms were expressed, only general inferences may prove useful. 
Terms from this type of equation will be used and tested in 
Chapter VI. 
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CHAPTER VI 


AVE-IV STABILITY ANALYSES 
A. Introduction 

This chapter will deal with the AVG and LAG AVE-IV profiles 
as they are applicable to the current standard stability indices of 
thunderstorms and severe weather reported in Chapter IV. The profile 
indices will be compared and the ability of each to forecast/measure 
severe convective weather will be determined. Temperature- 
dependency for each index will also be included, as well as the 
introduction of a new lag index. 

B. Temperature-Dependence 

It is generally desirable to know if a stability index 
changes with the changing temperature of an air column. The index 
is said to be temperature-dependent if this is the case. An index 
with a very large temperature-dependency is undesirable for use in 
representing the stability over a large geographical area in which 
differing air masses may reside. The changing air characteristics 
would affect the threshold value of an index, as reported by 
McPherson [29]. 

McPherson [29] suggested a method to determine the 
temperature-dependency of an index. His method is the approximate 
procedure used in the present investigation. A range of stability 
conditions for each index is determined first. The selection of two 
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hypothetical vertical temperature profiles, which represent (a) near 
normal stability and (b) less stable conditions, is accomplished by 
assuming moist-adiabatic and dry-adiabatic vertical temperature 
lapse rate conditions, respectively. These two conditions are 
chosen because they separate the area for conditional stability, and 
also give index values on either side of the threshold values 
(Chapter IV) which are representative of thunderstorm/severe thunder 
storm conditions as used in this study. For both adiabatic 
conditions, the temperature at 700 mb was assumed fixed at 0°C, and 
the dew-point depression at all levels was assumed to be 10°C. The 
resulting range for all indices is presented in Table 19. 

Index values were next obtained given five different cases 
of moist-adiabatic lapse rate conditions from 900 to 500 mb, with 
0^ values of 0, 6, 12, 18, and 24°C, assuming saturation at all 
levels. These index results are presented in Table 20, which now 
gives a relationship between index value and temperature changes. 

The range of index values obtained from the five different tempera- 
ture (0^) cases is presented in the next-to-last column of Table 20. 
The right-most column of Table 20 presents a percentage change of 
the range of each index with respect to the index's total range, as 
given in Table 19 for the five temperature categories. 

From these calculations it is shown that the Showalter and 
Bradbury indices have small or no temperature-dependence since they 
involve conservative moist-adiabatic procedures (0^) in their compu- 
tation. The Jefferson, Modif ied-Jefferson, 0*, Energy, and Modified 
Martin indices are only moderately (~1 to 10%) temperature-dependent 
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Table 19. Stability Index Range Determined by Moist 
and Dry Adiabatic Lapse Rates 


Index 

Normal 

Stability 

Moist- 

Adiabatic 

Lesser 

Stability 

Dry- 

Adiabatic 

Range 

of 

Values 

SWEAT 

8 

549 

541 

Vertical Totals 

26 

41 

15 

Cross Totals 

16 

31 

15 

Total Totals 

41 

72 

31 

CD 

m * 

-14 

3 

17 

Showalter 

7 

-10 

17 

Rackcl iff 

28 

44 

16 

Jefferson 

34 

55 

21 

Mod. Jefferson 

21 

42 

21 

Boyden 

2.5 

2.5 

0 a 

Bradbury 

3 

-7 

10 

K-Index 

14 

37 

23 

Energy 

2 

-4 

6 

Mod. Martin 

10 

-20 

30 


a In defining the 700-mb temperature as a 
constant here for both adiabatic processes, the 
Boyden index will result in a constant value. 
McPherson [29] concludes that the Boyden index is 
very temperature-dependent. 
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Table 20. Stability Index Dependency as a Function of Temperature Change 


Index 

CJ> 

O 

o 

6°C 

6 

w 

12°C 

18°C 

24°C 

Range 

of 

Values 

0^ Range 3 as 
Percent of 
Table 19 Range 

SWEAT 

-76.6 

-71.0 

-36.2 

44.2 

123.4 

200 

37% 

## 

Vertical Totals 

32.6 

30.7 

26.2 

23.9 

20.6 

12 

80% 

## 

Cross Totals 

22.6 

20.7 

16.2 

13.9 

10.6 

12 

80% 

## 

Total Totals 

55.2 

51.4 

42.4 

37.8 

31.2 

24 

77% 

## 

9 e 

2.8 

3.0 

2.7 

3.2 

3.6 

0.9 

5% 

# 

Showalter 

0 

0 

0 

0 

0 

0 

0% 

* 

Rackcliff 

41.4 

38.7 

33.3 

30.3 

26.4 

15 

94% 

## 

Jefferson 

41.4 

42.3 

40.5 

41.1 

40.8 

1.8 

9% 

# 

Mod. Jefferson 

28.4 

29.3 

27.5 

28.1 

27.8 

1.8 

9% 

# 

Boyden* 5 

22.7 

14.9 

6.9 

-1.0 

-8.6 

31 

— 

## 

Bradbury 

0 

0 

0 

0 

0 

0 

0% 

* 

K-Index 

3.8 

8.7 

11.1 

15.5 

18.8 

15 

65% 

## 

Energy 

1.21 

1.28 

1.80 

1.60 

1.55 

0.59 

10% 

# 

Mod. Martin 

-1.2 

-1.6 

-2.0 

-2.4 

-2.6 

1.4 

5% 

# 


Temperature-dependency code: * = None or small; # = Moderate; ## = High 


McPherson [29] concludes that the Boyden index is very temperature- 
dependent. 


IIIIIIIIII 


Finally, the SWEAT, K, three Totals, Boyden, and Rackcliff indices 
are all highly dependent upon temperature changes. Since 0^ was 
derived by using the 0g £ approximate equation for 0^, the resulting 
moderate temperature-dependency actually would have been less had 
this equation not been used. 

Temperature-dependency on an averaged profile, as compared 
to an individual provile, may appear unimportant since only averaged 
(AVG or LAG) severe storm thermodynamic profiles and indices are 
developed. However, when other individual vertical soundings (and 
their computed stability indices) are compared to those of the 
averaged profile, temperature-dependent indices may give unrealistic 
results. Also, the results of the AVE-IV soundings being averaged 
springtime soundings, over two independent days for the east/central 
United States, may not apply accurately for a different season 
(temperature regime) or location, if temperature-dependent indices 
are used. Therefore, in the evaluation of index performance, it is 
well to know in advance which indices are temperature-dependent and 
which are not. 


C. Stability Index Results 

This section presents the stability index results, based on 
the LAG and AVG profiles. The 14 stability indices described in 
Chapter IV were used in conjunction with the MDR LAG and AVG 
averaged atmospheric profiles. These results are presented in 
Table 21. 
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Table 21. Stability Index Values for LAG and AVG MDR Profiles 


Index 


MDR- 

-LAG 



MDR- 

-AVG 


Approximate 

Threshold 

Index 

Value 

A 

B 

C 

D 3 

A 

B 

C 

D a 

SWEAT 

194 

237 

271 

C331^> 

186 

233 

249 

290 

250 to 350 

Vertical Totals 

26.6 

26.1 

27.3 

29.4 

26.6 

26.0 

27.6 

(guT) 

26 

Cross Totals 

17.7 

20.9 

22.0 

CH3> 

17.1 

21.1 

22.0 

22.1 

18 

Total Totals 

44.3 

47.0 

49.3 

dUD 

43.7 

47.1 

49.6 

51.6 

44 to 50 

e*° 

^ r 

-3.5 

-3.2 

-6.0 


-2.9 

-3.4 

-5.3 

-8.1 

— 

Showalter 

4.1 

1.9 

0.0 

-1.2 

3.6 

1.7 

-0.6 

C2D 

-3 to -6 

Rackcliff 

29.1 

30.4 

31 .3 

<32J> 

28.8 

30.2 

31.3 

32.0 

30 to 35 

Jefferson 

38.1 

39.8 

41.7 

(JO) 

37.6 

39.5 

41.5 

42.9 

30 

Mod. Jefferson 

20.8 

26.2 

27.1 

( 29 T 0 ) 

19.6 

26.2 

27.8 

28.6 

28 to 29 

Boyden c 

6.8 

6.4 

5.8 

5.2 

6.7 

6.5 

5.2 

CUD 

— 

Q 

Bradbury 

0.8 

-0.1 

-2.0 

(2D 

0.7 

-0.4 

-1.5 

-2.4 

— 

K-Index 

11.5 

21.1 

22.2 

24.7 

9.6 

21.8 

23.9 

cgD 

30 to 35 

Energy 0 

0.20 

-0.27 

-1.08 

dD 

0.26 

-0.18 

-0.94 

-1.61 

<-2 

Mod. Martin 0 

2.2 

0.7 

-3.4 


2.6 

0.3 

-3.8 

-4.8 

— 


a Circled "D" category values indicate the largest unstable index value. 


b c 

Potential Lag index. Indices in which instability is negative (-). 


i limn 


As can be seen from Table 21, the category D profiles pro- 
duced the largest instability index values, as one would expect. 
Since category D conditions are the main items of interest in the 
present investigation, emphasis will be placed on them. The circled 
category D LAG and AVG stability indices indicate the index with the 
largest D-category index value. Two indices, in particular, show a 
much greater LAG instability than their AVG counterpart index value. 
The SWEAT and modif ied-Martin indices both indicate a LAG-REG 
difference greater than It, as a function of the index range. These 
two indices would be potential LAG stability index forecast indi- 
cators when used prior to the occurrence of severe weather. 

To establish concrete threshold values for all the indices 
used here is difficult, since each index may offer a threshold index 
value for only a selected type of thunderstorm condition ( i . e - , 
scattered thunderstorms rather than numerous severe thunderstorms). 
However, an attempt has been made to include approximate threshold 
index values for severe-type thunderstorms (see last column of 
Table 21). 

Note that most indices equal or exceed the threshold values 
indicated, with the exception of perhaps the K-index. However, the 
K-index has been designed for routine, non-severe thunderstorm pre- 
diction. 

Therefore, all indices presented in Table 21 appear to be 
potentially equal by this analytical comparison between prior and 
actual severe storm averaged conditions. This supposition will have 
to be considered in Chapter VII, when an actual, independent set of 


severe storm soundings is presented and analyzed with respect to 
atmospheric stability. 

D. Johnson Lag Index 

As contemplated in Chapter V, Section G, it is believed that 
the development of a forecast-type procedure or index should be 
attempted that is based entirely upon the differences noted in the 
averaged AVG and LAG profiles. If the environment three to six 
hours prior to severe weather shows any type of parametric structure 
difference from that at the time of severe weather, a stability 
index/procedure should be developed to model this phenomenon. Since 
wind differences are small between LAG and AVG profiles, and the 
individual wind profiles are so variable, it was felt that for the 
initial attempt, winds should not be used — only the use of signifi- 
cant thermodynamic parameter changes versus altitude to keep the 
index simple. 

As explained earlier, the major differences observed in the 
temperature structure between LAG and AVG profiles occur throughout 
the 800- to 650-mb and 650- to 500-mb levels. The main 0^ differ- 
ences noted occur between the 900- to 800-mb and 750- to 700-mb 
levels. The LAG and AVG temperature and equivalent potential 
temperature lapse rates that exist between these pressure levels 
were then calculated. A gradient halfway between the LAG and AVG 
gradients was selected as being a most representative standard of 
atmospheric conditions between three hours prior to storms and storm 
occurrence itself. Lapse rates on one side of this standard gradient 


99 


would represent conditions of the LAG, while gradients observed on 
the other side of this standard would represent AVG conditions. 

The four thermodynamic terms mentioned earlier were therefore 
selected as potential forecast terms: Two terms to represent 

temperature gradients in lower and upper atmospheric areas, and two 
0£ gradient terms to represent the low- and middle-atmosphere temper- 
ature and moisture structure. The four terms were then combined so 
as to maximize the negative value of the index in representing 
.extreme instability only during LAG-D time (three to six hours before 
storms). Since this gradient procedure, or index, is maximized a 
few hours before storm occurrence, the application of the index 
during periods of severe weather (AVG-D conditions) should result in 
a positive value. This new Johnson Lag Index (JLI) is expressed as: 

JLI = (-11.5 - a T 6 50-800^ + 2 ^ AT 500-650 + 14 ' 9 ^ 

+ 2 ( A0 E 800-900 + 3 ’ 5 ^ " 1/3 ( 3 * 0 + A0 e 700-750^ ’ ^ 

where. 


650-800 

= T 650 ' 

T 800 

9 

500-650 

= T 500 

T 650 

9 

'E 800-900 

0 E 800 

- 0 E 

900 ’ 

'E 700-750 

= 0 E 700 

“ 0 E 

750 ’ 


(T and 0£ units in °C or °K) . 
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The four terms of Eq. (25) were weighted by applying multi- 
plication factors of 1, 2, 2, and 1/3, respectively. This was done 
to offset the effect of the category A (non-precipitation) small 
temperature and potential temperature gradients, which tended to 
allow the unweighted JLI equation to produce an unstable negative 
JLI value close in magnitude to LAG-D JLI conditions. Thus, this 
weighting will help eliminate the occurrence of false alarms when- 
ever category A, non-precipitation areas are encountered. The 
weighting factors were determined from a subjective, tri al -and-error 
procedure involving different combinations of weighting, in order to 
arrive at a large JLI difference between A and D precipitation 
conditions. 

The JLI values calculated for LAG-D conditions equaled -4.35. 
Likewise, JLI values computed for AVG-D conditions resulted in a 
value of +2.76. The theory, then, is that if atmospheric conditions 
from an individual sounding produce a negative JLI of similar or 
greater magnitude, one should expect severe weather to occur within 
the next three to six hours. This conclusion has yet to be proven, 
and is only stated at this time. In Chapter VII the theory will be 
tested as to its performance along with the other stability indices. 
The complete JLI values versus MDR categories of LAG and AVG profiles 
are given in Table 22. 

Again, it should be realized that the very small parametric 
differences noted in these averaged profiles have been used in the 
construction of the JLI; whereas, in reality, individual atmospheric 
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Table 22. JLI versus MDR Categories for 
LAG and AVG Conditions 


MDR 

Category 

LAG 

JLI 

AVG 

JLI 

A 

0.52 

-0.15 

B 

4.11 

4.78 

C 

1.78 

3.45 

D 

-4.35 

2.76 
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soundings do, indeed, have a much greater range of variability in 
the vertical. The question is, how well will the JLI model the real 
atmosphere? 
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CHAPTER VII 


STABILITY CRITERIA APPLIED TO AVE-SESAME-I 
A. Introduction 

This chapter will present an analysis which will use the 
stability indices defined in the previous chapters. These are 
applied to a different and independent set of individual data 
soundings taken during severe weather situations. The 
AVE-SESAME-I [2] data case of April 10-11, 1979 was selected for the 
present study as the comparison data set against which to run all of 
the indices. 

The synoptic situation for AVE-SESAME-I will be presented 
along with the individual soundings. Stability indices will be 
computed for all soundings prior to, during, and after severe storm 
occurrence. The computed indices will be compared in helping 
determine how each index varies throughout this data set, and how 
each index might be used as a short-term predictor of severe weather. 

B. Synoptic Situation 

The AVE-SESAME-I time period was chosen for the stability 
index evaluation case because the AVE-IV and AVE-SESAME-I projects 
involved April storm cases in which similar synoptic weather situ- 
ations developed. 

A low-pressure system located north of western Texas with 
associated frontal positions existed, allowing a moist Gulf flow to 
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persist over the southeastern and southern/middle plains areas of 
the United States. This situation, coupled with the advancing cold 
front, caused extensive convective and severe weather to form, with 
the development of two pre-frontal squall lines, during both AVE 
cases. The destructive Wichita Falls, Texas, tornado that occurred 
during AVE-SESAME-I was only one of more than 40 tornado occurrences. 

Abilene, Texas, was the sounding station chosen to analyze 
during AVE-SESAME-I because it was the closest station during most 
of the tornado and severe weather occurrences in the north-Texas and 
southern-Okl ahoma areas. A preliminary weather summary for 
AVE-SESAME-I has been published by Williams [68], and the individual 
upper-air sounding data are available in the Gerhard et al . [2] 
document . 

Since the Abilene sounding data were selected for analysis 
in this study, the time and location of the severe weather /tornado 
occurrence around Abilene during April 10-11, 1979 will be discussed 
at this time. Generally speaking, there were three separate severe 
weather patterns which occurred near Abilene during the afternoon of 
April 10, 1979 and extended through the early morning hours of 
April 11, 1979. 

The first severe weather event consisted of hail damage 
between 125 and 150 miles north and northeast of Abilene, between 
1730-1800 GMT on April 10. The second very severe weather outbreak 
occurred between 2050-0100 GMT, with tornado and hail occurrence 
from 75 to 150 miles north and then north-northeast of Abilene. 

This included the Wichita Falls tornado. Finally, a squall line 
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developed around 0245 GMT (April 11) from Abilene and extended 
~75 miles south-southwest to San Angelo, Texas. For the next six 
hours (until 0817 GMT), this squall line moved eastward producing 
hail and some tornadoes from 35 to beyond 125 miles of Abilene. The 
severe weather occurred at Abilene and then moved south of, southeast 
of, east of, and finally, northeast of the city. Figure 27 indicates 
the severe weather pattern which occurred around Abilene (ABI) from 
1200 GMT on April 10 to 1200 GMT on April 11, 1979. 

C. Sounding Analyses 

For the Abilene site, eight atmospheric soundings were taken 
between 1200 GMT on April 10, 1979 and 1200 GMT on April 11, 1979. 
Only the 0600 GMT sounding of April 11 was missing due to tracking 
problems, with no second release as backup. 

Five of the eight critical Abilene severe weather soundings 
are presented in Skew-T form in Figures 28 and 29. Given are the 
April 10, 1979 soundings for 1442, 1740, 2034, and 2333 GMT, together 
with the one for 0226 GMT on April 11, 1979. The progression of 
these soundings in time indicates that low-level moisture was con- 
fined by a capping inversion to levels under 800 mb prior to 1442 GMT. 
From these profiles, extremely dry air can be seen above this level. 
After 1500 GMT, the capping inversion lifted and, with storm develop- 
ment, allowed moist air to penetrate upward to beyond 600 mb by 
2034 GMT. The April 11, 1979 Abilene sounding for 2333 GMT shows 
that during this time period, while severe weather was currently 
affecting the Wichita Falls area, the dry-line behind this system 
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Figure 27. Severe weather occurences between 1200 GMT on April 10, 
and 1200 GMT on April 11, 1979 in the south central 
United States [68]. 
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Figure 28. Abilene, Texas, severe weather soundings for April 10, 1979. 
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had now moved into the Abilene area. This condition was relatively 
short-lived because a second squall line was forming near Abilene 
and started inflicting severe weather there and eastward by 0245 GMT. 
The 0226 GMT sounding of Figure 29 shows an abrupt increase in 
moisture up to 350 mb, where the data terminate. Sometime after 
0600 GMT on April 11, 1979, the cold front began entering the 
Abilene area, bringing dry air close to the surface while still 
leaving a pocket of moisture above the front between 550 and 750 mb. 

Since stability is the item of interest in the present 
investigation, the 15 stability indices used earlier were computed 
for each Abilene sounding. These stability index results are pre- 
sented in Table 23, together with the exact time of radiosonde 
release. Listed below the index values in this table is a severe 
weather timeline applicable to the north-central Texas area, within 
150 miles of Abilene. This separation of sounding site and area of 
severe weather occurrence is, indeed, too large to be completely 
applicable to the Abilene data. Therefore, one should keep in mind 
that the Abilene timeline needs to be shortened somewhat. Also on 
Table 23, the highest three unstable index values for each index 
have been circled for easy reference. The most unstable value has 
also been marked with a superscript "a." 

As can be seen in Table 23, there seems to be good general 
agreement that most all indices appear to perform adequately in the 
evaluation of atmospheric instability during the passage of the two 
squall systems near Abilene. Profiles 4 and 6 were the two soundings 
taken at Abilene just prior to the severe weather which occurred 
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Table 23. Abilene, Texas, AVE-SESAME-I Sounding Stab i 1 ity Index Values 


April 10, 1979 


April 11, 1979 


Sounding No. 
Time (GMT) 
Index 


1121 


1442° 


1740° 


2034 


2333° 


0226° 


0600 


7 

0806 


8 

1105 


SWEAT 221 69 C472~) 

Vertical Totals 28.5 (50) 27.8 



Boyden 

Bradbury 0 

K-Index 

Energy 

Mod. Martin c 
JLI C 


-3.2 

-0.7 

0.2 

- 0.2 

-7.3 

^29 



Abilene Area Severe Weather 
T imel ine GMT : 


Description: 


(No Convective 
Activity) 


1730- 
1800 
(Hail ) 


2050 


0100 


0245 


-33 

26.6 


-7.1 

2.8 

15.3 

1.9 

13.5 

6 


0817 


(Tornadoes 
and Hail) 


(Hail and 
Tornadoes) 


-37 

26.3 


Cross Totals 

16.1 

3.1 

CO J> 


15.0 

dO) — 

9.3 

6.3 

Total Totals 

44.6 

36.0 

<53J> 

(sO) 

50.0 

(62j£ — 

35.9 

32.6 

e* c 

-11.3 

-8.1 


CJ6J> 

-6.0 

<S3> — 

11.2 

14.0 

E c 

Showal ter 

-1.1 

7.0 

(ED 


0.8 

— 

8.9 

10.3 

Rackcl iff 

32.2 

C33J) 

32.6 

(50) 

31.4 

— 

26.0 

26.1 

Jefferson 

43.0 

44.3 

C44T?) 

( 4 O) 

42.4 

(50) 

33.1 

32.6 

Mod. Jefferson 

19.9 

20.8 

10.5 

( 50 ) 

( 2 O) 

(fO) — 

21.9 

22.9 


-6.9 

3.8 
13.9 
2.3 

13.8 
12 


(Storms Move 
Eastward) 


Circled values indicate the highest three unstable index values for each index. 
b Most unstable stability index value. indices in which instability is negative (-). 




near and around the city. Stability index values from Table 23 
indicate that most indices peak (with instability) using soundings 4 
and 6 data; 10 of the 15 indices peak using sounding 6, while three 
peak using sounding 4. This means that 13 of the 15 peaked during 
the occurrence of upper-level moisture buildup, just prior to the 
onset of the Abilene storms. Only two indices (0^ and JLI) peaked 
at times prior to this. Sounding 6 is more unstable than sounding 4 
because the storms developed very close to the sounding site, and 
the moisture aloft had developed more extensively than during 
sounding 4. The dry-line passage at Abilene between 2200-0000 GMT 
can readily be seen by the sudden increase in stability in most all 
of the indices during sounding 5 (2333 GMT). While weather activity 
existed eastward of Abilene during sounding 8 (0806 GMT, April 11), 
all indices show a general increase in stability as the cold front 
arrives. 

Table 23 also hints that soundings taken when storms are not 
in progress in the general area result in slightly greater insta- 
bility than when storms have formed in the area during the radiosonde 
release. This may seem to indicate that the instablity (stored 
potential energy) which can build up prior to storm occurrence can 
be relieved (made more stable) through the release of thunderstorm 
kinetic energy activity. 
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D. Exceptions to the Norm 


There are a few stability indices which peak at an earlier 
sounding than the rest. It was decided to look at each index that 
fell into this category. 

Since the 0* index peaked out during sounding 3, the cause 
was sought. The index is based on a 0^ difference between 700 mb 
and the surface level. Since sounding 3 indicated that the atmos- 
phere dried out very quickly between 750 and 700 mb, this would also 
produce a very dramatic 0^ drop between the same two levels, 
resulting in a very unstable 0* index value. 

The modified Martin index indicated a slight instability 
peak using sounding 2. This index uses the 850-mb level as the 
comparison altitude. At this time interval there existed a large 
capping temperature inversion top located at 850 mb. This large 
temperature value would produce a higher index value than if this 
cap inversion top were located at a different level. 

The Energy index peaking during soundings 3 and 4 is 
believed to be due to the ample abundance of moisture at 850 mb 
during these two sampling times, which would strengthen the 850-mb 
E-j- value. Sine the JLI was designed to peak out during time periods 
prior to storm development, this early peaking of the JLI index 
during soundings 1 and 2 is expected and desirable. 
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E. Lag Testing 


In Chapter VI it was indicated that, based on the AVE-IV LAG 
profile as it related to the AVG profile, three indices appeared to 
be potential lag indices. These three indices were: SWEAT, modified 

Martin, and the JLI. According to the AVE-SESAME-I sounding data 
(Table 23, page 111), all of these indices, with the exception of 
the JLI, fail to qualify as a lag index, since the peak index out- 
liers which occur before storm development have been explained away. 

The JLI does give large negative values (-29 and -18) during 
the non-storm time period represented by soundings 1 and 2. When 
distant storms occur, Abilene sounding 3 records a JLI = -7. Just 
prior to the first major outbreak of storms closer to Abilene, 
sounding 4 gives a JLI = -10. The dry-line passage, during 
sounding 5, produces a JLI = +1. Sounding 6, released 19 minutes 
prior to hail occurrence near Abilene (51 minutes prior to first 
tornado report) gave a JLI = -28. This large negative index value 
was surprising, since the sounding represents squall line-produced 
activity. However, the JLI could still be sensing the intense, 
unstable, pre-squall line environment which appears not to have 
passed the release site at this time. Overall, the JLI has 
functioned well and it gives large positive values (+6 and +12) when 
the cold front moved into the area. This indicates that no more 
storms were due to follow. 
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Based on only one severe storm case, it appears that of 15 
stability indices tested as a pre-storm lag index, only the JLI 
appears to give satisfactory results thus far. However, since the 
JLI is a new index, representing low- and middle-level temperature 
and moisture, it will have to be tested further, and possibly be 
adjusted, before it can quality as a lag/forecast index for severe 
storms . 
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CHAPTER VIII 


CONCLUSIONS 

The following project goals and conclusions were accomplished 
and presented in this study: 

1. Averaged AVE-IV environmental thermodynamic/wind profiles 
have been presented for different MDR severe storm 
periods (AVG), and for periods three to six hours prior 
to severe storm occurrence (LAG). 

2 . The AVE-IV AVG and LAG profiles were analyzed paramet- 
rically and with 14 common atmospheric stability indices 
to determine if a severe storm forecast index or pro- 
cedure could be developed based on these averaged severe 
storm profiles. 

3. A thermodynamic lag index, called the Johnson Lag 
Index (JLI), was developed, based upon low- and middle- 
tropospheric level temperature and moisture structure 
using the AVE-IV averaged profiles. The JLI was designed 
to have short-termed forecasting ability. 

4. Based on the averaged AVE-IV profiles, two other 
stability indices (SWEAT and modified Martin) had some 
potential as forecast lag indices. 

5. All 14 stability indices and the Johnson Lag Index were 
tested by employing an independent severe storm case 
study using the AVE-SESAME-I individual data soundings 
from one station. 
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6. 


All AVE-SESAME-I stability indices tested appeared to 
recognize the severe weather environment with unstable 
values, as well as presenting stable values when severe 
weather had passed. 

7. All AVE-SESAME-I stability indices tested as a pre-storm, 
three-hour lag forecast index performed unsatis- 
factorily. Only the JLI appeared to show promise in 
terms of forecasting severe weather three to six hours 
prior to occurrence. However, more testing of this 
index with case study data is needed. 
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